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PREFACE 

This report is the result of the "TYCHO" Study Group Summer 
Study Session held at Dartmouth College, Hanover, New Hampshire, 
during the period June 20 to July 21, 1967. 

The report is organized into a basic section and appendices 
containing a series of scientific papers written by individual 
"TYCHO" Study Group Members as research contributions. The basic 
section contains recommendations, partially supported by the 
appendices, concerning NASA’s proposed planetary explorations 
program. 

The 196 7 "TYCHO." Summer Study Session investigations were 
confined to matters relating to atmospheric and surface conditions 
of the planets Venus and Mars of importance to proposed manned or 
unmanned landings . The investigations consisted of a comprehen- 
sive analysis of published data and presentations by the best 
available experts in the field. 
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REPORT 


OF 

1967 "TYCHO" SUMMER STUDY SESSION 
June/ July 1967 
SUMMARY AND RECOMMENDATIONS 

During July, 1967, the "TYCHO" Study Group met at Dartmouth 
College to examine various aspects of the NASA program for 
planetary exploration, Considerations of the group indicated 
that in the near future greatest scientific interest would be 
associated with unravelling the nature of the atmosphere of 
Venus and the details of the surface of Mars. The major impetus 
for planetary exploration stems in large measure from the 
opportunities to shed light on the origin of the solar system 
and the origin of life. These two questions were considered by 
the National Academy of Science and discussed in the Woods Hole 
Report. In response to this report, NASA proposed the planetary 
exploration program shown in Fig. 1. This program amply responds 
to all of the issues raised. Recognizing the completeness of 
both the report and the outlined program, the present study group 
raised the question, "What should result if delays arise as a 
result of either fiscal or technical problems?" In response to 
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Figure i Evolution and Development of Planetary Missions. 

Following the Voyager missions to Mars, starting in 1973, will be 
the Voyager missions to Venus, starting in 1977. In the 1980's, 
Voyagers will be used to explore the outer planets. It is important 
to note that, during the period in which the large Voyager system is 
first used for Mars, there is a continuing need for Mariner-class 
spacecraft for precursor missions to Venus and to the outer planets. 
Within the restrictions of launch opportunities and the lead times 
needed to develop and test spacecraft, the sequence of missions 
shown remains highly flexible. Scientific and technical information 
gained from previous flights feed into subsequent missions. (From 
"Summary of the Voyager Program," NASA unnumbered publication, 
January 1967.) 
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this question, the following recommendations are made: 

Recommendat ion 1 . -F ollowing the successful Mariner fly-by 
mission (1969) and the Mariner atmospheric probe mission (1971) 
the 1973 opposition of Mars should be used (as planned) for an 
orbiting-landing vehicle (Voyager). In the event that the 
Voyager 1973 spacecraft will not be available, an additional 
Mariner spacecraft of the 1969 type should be used in an 
orbiting mode . This mission should be designed to supply 
multiple occultation experiments, high resolution optical sur- 
face mapping with a TV camera system, and an infrared mapping 
of the surface. The data obtained could provide information 
on the seasonal and diurnal fluctuations and give a better 
understanding of the atmospheric effects than would be possible 
with a second atmospheric probe. Furthermore, the multiple 
passes would allow the radio occultation experiment to provide 
some information on the contour of the terrain. 

Recommendation 2 . -R ecognizing the great uncertainty associated 
with the density and composition of the atmosphere of Venus, 
the next mission (1972) to Venus should be (as is planned) an 
atmospheric probe. In the event of delays or shortages in the 
Mariner program, the Mariner spacecraft/atmospheric probe 
scheduled for Mars in 1971 should be rescheduled for the 1972 
Venus opposition . 

Recommendation 3 . -Recognizing that the contributions from 
ground-based astronomy and in particular radar astronomy have 
been of significant importance to the space program (e.g. , 
the definition of the astronomical unit, orbit elements of 
the planets , rotation rates of Venus and Mercury, terrain 
elevation of Mars, and mean surface slopes of the planets), 

NASA should take positive steps to ensure the continuation of 
those radar activities now being carried on as peripheral 
activities at installations supported by the Department of 
Defense and other agencies. 
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The recommendation concerning an orbiting spacecraft about the 
planet Mars reflects the view that the 1964, 1969, and 1971 
Mariners will give enough information required for successful 
landing on the surface. However, in order to put a lander to best 
use, detailed information concerning the character of the surface 
will be needed. The 1973 orbiting spacecraft with a TV camera 
system would yield this information and allow examination of many 
possible landing sites. Additional studies of seasonal fluctuations 
of the dark and light areas would also be possible. (For detailed 
recommendations concerning the Mars' probe, see Appendix B.) 

The principal question related to Venus is the constitution 
and pressure of the lower atmosphere and its interaction with the 
surface. The current (1967) Venus fly-by can be expected to yield 
some information concerning the atmospheric density. However, 
until the composition of the clouds and the surface pressure is 
determined our knowledge of conditions there will remain extremely 
incomplete, and it will be impossible to make models for the 
physics and- chemistry of this planet. For this reason the 1972 
mission to -Venus should include an atmospheric probe. In light 
of the success of the 1964 Mariner IV occultation results and the 
fortunate coincidence that nearly the entire atmosphere can be 
attributed to one component (CC^), the need for an atmospheric 
probe of Mars is not pressing. Thus, it may even be desirable to 
reschedule the 1971 Mars /Mariner spacecraft probe to 1972 Venus. 

(For details concerning the Venus atmospheric probe, see 
Appendix G . ) 

Various particular questions are considered in Appendices 
A through J. 
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DISCUSSION 


The impetus for planetary exploration stems in large measure 
from the opportunities this program offers to shed light on the 
origin of the solar system and the origin of life. In the pro- 
cess of answering these questions a necessary first step is 
determination of the physical environment. At the present time, 
two most important scientific questions are the possible presence 
of surface or subsurface water on Mars and the gross features of 
the composition and the pressure of the atmosphere of Venus. 


Mars 

In recent years the study of Mars by optical and infrared 
spectroscopy, radar reflectivity, and the Mariner IV space probe 
has added much to the knowledge of observational astronomy, but 
little to a detailed understanding of the geology of the surface 
or interior. Factors now considered well determined are the 
atmospheric pressure, gross composition of the atmosphere, and 
the surface temperature. (See Appendix F). Of those Martian 
features apparent even to the early optical observers — dark and 
light areas, polar caps and dust storms to say nothing of the 
fabled canals — virtually no information has been added within the 
last decade. The Mariner IV mission radio occultation experiments 
successfully determined the number density of the atmospheric 
particles as a function of altitude. The same Mariner IV pictures, 
however, added an additional unresolved question. This result is 
the occurrence of large craters with an absence of sharp features. 
The implied erosion seems too effective to have been caused by 
winds alone. 

The question of the presence of water on Mars is of great 
interest because it has an important bearing upon the problem of 
extra-terrestrial biology and because it is closely related to our 
understanding of the degree and mechanism of erosion on that planet 
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The range of pressures and of surface temperatures on Mars 
straddles the triple point of water: 6.11 mb and 273°K, and thus, 

in principle , one could expect to encounter ice , liquid water and 
water vapor in amounts varying with altitude, lattitude , time 
of day and year. While there is evidence for water vapor in the 
atmosphere there is no comparable information about ice or liquid 
water. Polar caps are probably solid CO 2 but permafrost or ice 
may exist at preferred locations under a rather shallow coverage 
of soil. Intergranular moisture and a thin liquid top layer of 
permafrost could provide a suitable millieu for bacterial life . 

The dark areas on Mars are known to be warmer than the bright 
areas and if the recent radar studies are correct in indicating 
that these areas are topographic depressions rather than high 
plateaus , then it is natural to assume that they have more sub- 
surface liquid water or ice than the rest of the planet. This 
difference may account not only for their higher radar reflectivity, 
but also for the darker color either following the biological 
hypothesis or according to some of the inorganic models. It is 
recommended that radar reflectivity and surface temperature of 
the dark areas be measured as a function of diurnal and seasonal 
variations of insolation. Such observations could be best made 
from a Martian orbiter but ground-based studies could also be 
important . 

The present knowledge concerning Mars allows a reasonable 
coherent picture to be constructed of the pressure temperature 
composition of the atmosphere. This model has been used to pre- 
dict global wind patterns. The next series of questions concerning 
the Mars planet are related to the properties and character of the 
surface : 

1) What mechanism is responsible for the erosion? 

2) What is the origin of the season of color changes? 
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3) What is the nature of the polar cap? 

4) What is the relative terrain of the dark and light areas? 

The above questions will not be answered until samples of the 
surface are available along with detailed studies of the topology. 
The implication of this statement is that a landing on the surface 
of Mars with a modest automated laboratory will ultimately be 
necessary. This laboratory should be self contained, able to 
maintain observations during seasonal changes and have the 
possibility for executing a simple chemical analysis as well as 
monitoring the local weather. The present ignorance of the sur- 
face terrain features suggests that prior to landing it will be 
necessary to obtain relatively high definition optical photo- 
graphs of a certain fraction of the surface. The fraction of the 
surface subjected to such examination should include areas subject 
to seasonal changes. It is because of the selection of a suitable 
landing site for the automated laboratory that the lander probe 
of an orbiting mission was recommended in the previous section. 


Venus 

The planet Venus presents a more interesting challenge to the 
scientific community of 1967 than does Mars. This difference 
arises not from the amount of knowledge concerning Venus but rather 
the lack of it. Optical observations of Venus reveal only a 
generally mottled disk usually interpreted as clouds. The origin 
and composition of the clouds — a topic of much debate — has not yet 
been resolved. Present spectrographic information has indeed 
revealed the presence of C0 2 and H 2 0. These observations and the 
infrared and radar temperatures have been used to imply the 
presence of ice crystals or water vapor as the dominant mechanism 
producing clouds. However, alternate approaches involving high 
surface pressure and perpetual dust storms can be offered to present 
a consistent picture of the observed data. The conclusion from 
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these observations is that evidence obtained through spectroscopy, 
radar reflection, radio emission, or infrared temperatures will 
not yield a satisfactory model for the planetary atmosphere or its 
surface until the composition of the atmosphere is available. For 
this reason it is essential to program an early Venus fly-by 
equipped with an atmospheric probe which will determine the 
composition of the clouds and identify a composition profile of 
the atmosphere. 

The near equality of the radii, masses and uncompressed 
density of the Earth and Venus suggests that the composition and 
internal evolution of these two planets should be similar. The 
atmosphere of Venus is apparently quite different than that of the 
Earth even when the effects of geological activity of Earth are 
taken into account . Although still an unsettled problem the amount 
of H 2 0 in the atmosphere of Venus is apparently a great deal less 
than exists on the Earth. If the atmospheres of the Earth and 
Venus are secondary, i.e. due to outgassing of the interior, we 
would expect the amount of water outgassed from Venus to be simi- 
lar to the amount of water outgassed from the Earth if the 
primitive compositions of the planets were the same. 

The lack of water on Venus, if verified, could be due to: 

1) Less water available in the primordial accreting material 
from which Venue accreted as H 2 0 or as water of hydration in 
silicates. This is a possibility since Venus presumably 
accreted nearer the Sun than the Earth and the primitive 
particles were, therefore, slightly hotter. 

2) Less hydrogen available in Venus (relative to carbon) for 
reduction purposes. 

3) More effective removal of H 2 0 from the atmosphere of 
Venus by photodissociation and loss of hydrogen. 

If the latter possibility is important, then the lack of 
water in the atmosphere of Venus becomes less cosmologically 
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important. One can ask, what would happen to the Earth’s 
atmosphere if the Earth: 

1) Were moved to the orbit of Venus 

2) Had its axis of rotation less tilted to its plane or 

revolution, thereby decreasing seasonal effects 

3) Had its period of rotation increased to 243 days 

4) Had its magnetic field removed 

Would the cloud cover increase, the surface temperature rise 
to 600°K, and would coldtrapping effects be less effective in 
retaining water? These are essentially the observed differences 
between the Earth and Venus, and it is important to decide whether 
these differences are meteorological rather than cosmological. 

Some of these questions can be tackled by methods of 
theoretical meteorology, but a Venus atmosphere probe would tell 
us immediately the composition of the Venusian atmosphere and would 
allow us to compute the rate of escape of these components or their 
dissociation products. 

This is just one aspect of the problem of the origin and 
evolution of the solar system. The differences in composition of 
the Earth, Venus, Mars, Moon, Mercury, and the meteorites must be 
explained in any theory of the solar system, and it is important 
that differences which arise during evolution of the planets, such 
as the possibility of increased ^0 loss from Venus, be separated 
from primordial differences in the composition of the protoplanetary 
clouds . 
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ON THE RADIO OCCULTATION METHOD 
FOR STUDYING PLANETARY ATMOSPHERES 

ABSTRACT 

The problem of determining the refract ivity profile of a 
planetary atmosphere from optical or radio occultation data is 
identical in principle to the problem of determining the varia- 
tion of seismic velocities in the earth from the travel times of 
seismic body waves observed at the surface of the earth. In 
either case a complete set of data can be inverted uniquely, the 
only constraints being those fundamental to geometric optics. 

Expressions are given for converting radar doppler data to 
the index of refraction as a function of depth in the atmosphere, 
and for the effects of geometric spreading and absorption on the 
amplitude of the signal. 
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ON THE RADIO OCCULTATION METHOD FOR 
STUDYING PLANETARY ATMOSPHERES 
D, L. Anderson and R. A. Phinney 

July, 1967 
INTRODUCTION 

The purpose of this note is to point out a direct method of 
transforming radar doppler data from an occultation experiment 
directly into a profile of refractive index versus radius. 

Fjeldbo and Eshleman (1965) discuss in some detail the types of 
information contained in the phase and amplitude of a radio signal 
occulted by a planet. Their method of determining the refractive 
index profile is approximate but it should give good results 
without iteration for a planet with a thin atmosphere. They 
applied their technique to the interpretation of the Mariner IV 
occultation experiment on Mars. Kliore, et al. (1966), Kliore , 
et al. (1965a), Cain, et al. (1965), Kliore , et al . (1965b), 
Fjeldbo, et al. (1965), Fjeldbo, et al. (1966a, 1966b). For a 
thick atmosphere where strong bending of the rays is expected, a 
more general method of interpretation is desirable. Such a 
method, due to Herglotz, Wiechert, and others, has been well 
developed in the seismological literature. 

The Herglotz-Wiechert method is used to determine velocity- 
depth profiles in the earth from the observed variation of travel 
time with arc distance between source and seismometer. We give 
here a brief derivation in a form appropriate to the radio occul- 
tation experiment and a discussion of its applicability. For 
further discussion of such effects as focusing and defocusing, 
shadow formation, etc., we recommend Bullen’ s (1963) text. 

The geometry of the occultation experiment is summarized in 
Fig. 1. The coordinates are referred to the center of mass of 




the planet. The refractive index, n, is a function of ' r in a 
shell R < r < r Q where r Q is an arbitrary level lying effectively 
above all the mass of the atmosphere, and R is the radius of the 
solid planet. The ray path is straight-line outside the atmos- 
phere. It is assumed that the trajectories of the spacecraft and 
the earth are known with respect to the center of mass of the 
planet. The indicated spacecraft trajectory is the projection 
of the actual trajectory in the plane defined by the centers of 
mass of the spacecraft, earth, and planet. If the properties of 
the atmosphere depend only on r, then only components of velocity 
in the plane of the figure contribute to the observed Doppler 
shift . 

This research was performed as part of the 1967 "TYCHO" 
Summer Study at Dartmouth College under NASA Contract No. 
NSR-24-005-047 with the University of Minnesota. The authors 
are indebted to G. Fjeldbo for his discussion of the problem 
with the "TYCHO" group. 


REDUCTION OF THE DATA 


The observed Doppler shift is corrected for the relative 
motion of the earth-based receiver and the planetary center of 
mass. The residual Doppler shift is then converted into a total 
phase by 

4>(t) = 2 7rt Af ( t ) (1) 

where the zero of time is taken shortly prior to occultation. In 
the ray-optical description of the experiment, one can set up a 
family of constant phase surfaces and orthogonal rays. The 
selection of a time zero provides a unique value to the phase 
function as a function of the spatial coordinates. The space- 
craft trajectory provides a sampling of the phase along an arc. 
This information determines the angle of emission of the ray from 
the spacecraft and, eventually, the refractive index profile. 
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The eiconal equation for the phase is 



( 2 ) 


Assuming the spacecraft is above the atmosphere 

| | = w/c (3) 

From the Doppler data, the directional derivative — |- is 
determined along the spacecraft trajectory. The angle between the 
ray and the trajectory then follows: 


cos (a+t) = ||/|vt| * If (■*) 

where f may without error be taken as constant during the 
occultation. Since $ is known, a follows immediately. 

The ray is now projected to a reference sphere r . The 
angle of incidence on this surface is : 

sin i = r sin a/r (5) 

o s o 

(the subscript s refers to spacecraft coordinates). Solution of 
the trianges in Fig. 1 gives the entry point angle 


sin (£2 -£2 ) = [r cos a-r cos l ] (6) 

o s r s o o 

o 

and the angle subtended by the portion of the ray in the 
atmosphere: 

© = TT-i'-£2 (7) 

o o 

The rays to earth are sufficiently parallel that we can take 
• ^ • 

1 o ” o* 

Instead of Doppler as a function of spacecraft position, we 
now have the angle of incidence at the top of the atmosphere as 
a function of 0. We now transform to a relation between refrac- 
tive index and radius, using the method described by Herglotz 
(1907), Bateman (1910), and Wiechert and Geiger (1910). 
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TRANSFORMATION TO A REFRACTIVE INDEX PROFILE 


A ray parameter, p, is defined 

p = nr sin i (8) 

which is constant along a ray, by Snell's law, and serves as a 
label of the ray. In the absence of an atmosphere p would be the 
radius of closest approach, and may be called the impact para- 
meter of the ray. We have immediately: 

p = r^ sin i = r n (9) 

o o p p 

where the subscript p refers to the turning point of the ray. 

As an auxiliary variable, define ri , which varies along a ray, 

taking the value r at r = r and the value p at r = r : 

& o o P 

H = nr = p/sin i (10) 

Integrating along the ray gives the total arc subtended: 


0(p) 


2 I (p/r)(n 2 -p 2 ) 1/2 dr 
J r 

P 



, , w 2 2 , - 1/2 

(p/r)(ri -p ) 


dr 

dn 


dn 


( 11 ) 


The path length in the atmosphere is 

s = 2 J rn (n 2 -p 2 )" 1/2 dr 

P 

Now set 

2 2 
y = n -r Q 


w 


2 

P -r. 


2 


( 12 ) 


(13) 
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Then 

o 

0(p) = 2p J (y-w) -1/2 d ^ - dy (14) 

w 

This integral equation for In r(y) is in the form of the Abel 
integral equation and may be solved by reference to the Abel 
transform pair 

w 

f(w) = K J g' (y)(y-w)"^ 2 dy (15) 

o 

y 

g(y) = (Kir)”^ J f Cw) (w-y)”^ 2 dw 
o 


By inspection: 

y 

In [r(y)/r Q ] = ^ J 0 (w) 2p(w) "^(w-y) - ^^ 2 dw (16) 

o 

or 

n 

In [r(n)/r> ] = - f 

° tt J 



The right hand side is solved by the substitution p = n cosh q 
and integration by parts: 



The first term on the right hand side vanishes and exponentials 
can be taken: 
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r(n) = 


r Q exp { 


-* i 


cosh ^ (p/n ) ~ dp} 


(18) 


The integral is always positive; r Q > n and ^ is normally 
negative. It is possible to make ^ positive over certain 
ranges of p due to sharp gradients in n (n < 1), but the integral 
must always be positive. 

The refractive index profile r(n) follows from Eq. (18) and 
from the definition of n > Eq. (10). 


DISCUSSION 


Some of the hazards in applying Eq. (18) are discussed, in 
Jeffreys (1959) and Bullen (1963). Being due to the nature of 
the optical problem rather than the method of analysis, these 
considerations are equally applicable whatever method is used 
to derive the refractive index profile. In summary: 

1) d 0/dp may be positive for a limited range of p (Fig. 2). 
The function p(0) is then multivalued for certain ranges, due to 
multipath propagation. The spacecraft must be able to distinguish 
all multiple rays in order for an interpretation to be correct, 

2) In a neutral atmosphere, the increase of refractive index 
with depth causes the rays to be curved downward, In order that a 
ray emerge from the atmosphere by a direct refraction path the 
radius of curvature of the ray must exceed that of the tangent 
sphere. Since the downward curvature of a ray is given by 


d In n 

^ dr 

this condition is 

dn . n 1 

- ■— — < — ^5 — 

dr r ^ r 


(19) 


( 20 ) 
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P 


r o 

Figure 2 Schematic 0(p) Curves for Different Atmospheres. 

(a) No atmosphere: 0= 2cos - ^ (£— ) ; p > R. 

o 

(b) Strongly refracting neutral atmosphere 

(c) Weakly refracting ionosphere and weakly 
refracting atmosphere 

(d) Strongly refracting ionosphere and strongly 
refracting atmosphere 

(a) and (c) end occultation; (fc>) and (d) end asympototically 
at an optical shadow. 
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In the lower atmosphere of a planet like Venus Eq. (20) may not 
be satisfied, and rays with p below a critical value can emerge 
from the atmosphere only by reflection from the solid surface. 

The essence of this, or any quivalent inversion method is that the 
turning points (r^) of the rays considered must range continuously 
through the atmosphere. Since this is true only from r Q down to 
the critical level, Eq. (18) cannot be applied to the atmosphere 
below the critical level. In seismology, refractive index 
decreases at great depth result in the eventual emergence of the 
downward refracted rays, but this is extremely unlikely in 
atmospheric studies. 

A slightly different condition defines when a ray having the 
accessible range of p, 0 < p < r , can have a turning point at 
depth r: 

n r < r (21) 

o 

If a screening region of low index (ionosphere) lies above the 
neutral atmosphere, Eq. 21 is replaced by the more stringent 
condition. 


n r < n* r* 


( 22 ) 


where * refers to the ionospheric level of maximum nr. 

For an ionosphere with the scale determined for Mars by Mariner IV, 
8 3 

about 6 x 10 electrons/cm would be required. Consequently, 

Eq. 22 is unlikely to apply in a real case. 


The intensity variation of the refracted signal due to 
focusing and defocusing can be stated in the notation of this 
paper. Neglecting diffraction effects at caustics. 


I (exit) _ sin ifj 

d 0/ dp 

I(entry) sin(©+^) 

Ftan i c - tan a) sin a > tan -1 

V r Q cos ip r Q sin i Q 


( 23 ) 
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where \J> = S2 q - Q g , If r g is nearly equal to r , the expression 

in braces is small and positive; I (exit) will be very large 

(focusing) if^-=0 + {}«0. Ifr -* 00 , a -»■ 0 , -*■ i , and 

a P S ° 

focusing will occur if ^ = tan i Q . In the study of atmospheres, 

i Q is always near it/ 2 and that condition will not be satisfied 

(Fig. 2). 


Rays which are refracted downward at nearly the curvature 
of the planet have large ^ and are severely defocused. 
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VARIATIONS IN THE RADAR CROSS SECTION OF VENUS 


ABSTRACT 

Variations in the radar cross section of Venus have been 
reported by a number of observers. At decimeter wavelengths 
these are of the order of +_ 50% of the mean and have been 
attributed to variations in the nature of the terrain occupying 
the subradar point. At 3.8 cm the cross section is most 
frequently observed to be 'v .012% of the projected area of the disk, 
but values three times larger have been reported. The suggestion 
has been made that these increases result from changes in the 
amount of attenuation caused by the atmosphere of Venus. This 
would imply large scale variations in the distribution of the 
microwave absorbing agent and would have serious consequences 
for many of the atmospheric models that have been proposed. 

This paper reviews the evidence for variations in the cross 
section and provides revised values of the cross section at 
3.8 cm. The mean value for the cross section observed in 1966 
is now placed at .017. A closer examination of one instance in 
which the cross section at 3.8 cm was seen to increase reveals 
that this event was associated with the passage of a large 
feature through the region of the subradar point. This suggests 
that the variations observed at this wavelength, like those at 
longer wavelengths are associated with the nature of the terrain 
at the subradar point and are not the result of weather-like 
phenomena on Venus. 
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VARIATIONS IN THE RADAR CROSS SECTION OF VENUS 

J . V . Evans 
Lincoln Laboratory+ 

Massachusetts Institute of Technology 

INTRODUCTION 

In a companion paper we have reviewed the results of radar 
studies of the planets Venus and Mars. In that paper we 
discussed the variability in the radar cross section of Mars, 
but did not report on similar variations observed for Venus. In 
this report we review the evidence for variations in the cross 
section of Venus, and present hitherto unpublished data which 
throws light on the behavior at 3.8 cm wavelength. 

We are indebted to the staff of the Haystack Microwave 
Facility who were responsible for making the 3.8 cm measurements, 
to Dr. G. H. Pettengill for helpful advise and discussions, and 
to Dr. I. I. Shapiro and his colleagues for furnishing a 
tabulation of the longitude and latitude of the subradar point 
on Venus as a function of date. This work was performed during 
the course of the "TYCHO" Study Group meeting at Dartmouth College, 
Hanover, New Hampshire, June-July , 1967. 

METERWAVE AND DECIMETER OBSERVATIONS 

Measurements at the radar cross section of Venus have been 
made at meter wavelengths by James and Ingalls (1964) and 
Klemperer et al. (1964) at the wavelengths listed in Table 1. 

Both groups found that the mean cross section was a little under 
0.2 of the projected area of the disk, but that values as large 
as 0.6 were occasionally observed. 


+Supported by the U. S. Air Force 



TABLE 1 


CROSS SECTION OBSERVATIONS OF VENUS 


Observer Wavelength Period 


James and Ingalls (1964) 7. 
Klemperer et al. (1964) 6. 
Pettengill et al . (1967) 70. 
Evans et al. (1965) 23. 
Goldstein (1964) 12. 
Carpenter (1966) 12. 
Evans et al. (1966) 3. 


84 

iru 

11/6/62 

- 12/7/62 

00 

m. 

11/28/62 

- 12/7/62 

00 

cm. 

2/22/64 

- 9/20/64 

0 0 

cm. 

3/5/64 

- 8/26/64 

60 

cm. 

9/7/62 

- 12/15/6 

60 


5/10/64 

- 7/18/64 

80 


1/18/66 

- 6/27/66 


In the case of the A = 7.84 m measurements plane polarized 
signals were transmitted and received, and hence some of the 
observed variability must have been caused by Faraday rotation 
in the earth's ionosphere. An additional fading mechanism can 
be anticipated; namely, the constructive and destructive 
interference of signals reflected from different portions of 
the surface of Venus. Klemperer et al. (1964) show that the 
echo amplitudes have a Rayleigh distribution as would be expected 
due to this effect. Similar fading is encountered for the moon 
as its aspect changes with respect to a terrestrial observer, 
and has been studied extensively (Evans 1965). In the case of 
Venus near inferior conjunction the fading period would be of 
the order of _> 20 sec at A = 6 m so that in a five-minute 
observing period there would be 15 or less independent echo 
samples. Hence, much of the variability observed by James and 
Ingalls (1964) can be attributed simply' to the small number of 
independent signal samples obtained per run. Long period 
variability in the cross section may also have been present in their 
result, but is probably masked by these short period fluctuations. 
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Measurements of the cross section at decimeter wavelengths 
have been reported by Pettengill et al. (1967), Evans et al. 

(1965), Goldstein (1964), and Carpenter (1966) for the time 
periods and wavelengths listed in Table 1. At these wavelengths 
the fading period is reduced to the order of one second so that 
the statistical sample obtained during a single 5-minute 
observation period is considerably improved. In addition, all 
these observers employed antennas which could follow Venus, and 
thus were able to make repeated observations during the course 
of the day. In these measurements, therefore, the principal 
sources of uncertainty in the cross section were those introduced 
by poor signal-to-noise ratio and the calibration of the radar 
equipment . 

Goldstein (1964) obtained a mean value for the radar cross 
2 

section of .095 nr (where r is the radius of Venus) with 

departures from the mean as large as + 50%, which he attributed 

to variability in the scattering behavior of the planet. 

2 

Carpenter (1966) reported a mean value of .114 nr at the same 
wavelength, and found variations in this of the order of _+ 2 0%. 

No explanation for the difference between this result and 
Goldstein's has been offered. 

Evans et al. (1965) and Pettengill et al. (1967) found some- 
what larger variations at X = 23 and 70 cm respectively. Since 

the mean of the values published by Evans et al. is about 0.13 

2 . . 

nr , and the mean of the values obtained by Pettengill et al. 

2 

(1967) is .14 nr , we have compared the two sets of data m a 
single plot (Fig. 1) showing cross section vs date. Much of the 
variation shown in Fig. 1 must be introduced by the errors of 
measurement. These would be least in June and July (i.e. around 
inferior conjunction) and become progressively larger away from this 
time. Despite this the occasional values that are either half or 
twice the mean probably represent real changes in the cross 
section. 
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VENUS CROSS SECTION OBSERVED IN 1964 



Figure 1 Millstone Hill (X = 23 cm) and Arecibo Ionospheric Observatory (X 
cross section measurements as a function of date in 1964. 





The data shown in Fig. 1 span a time interval longer than the 
axial rotation period of Venus with respect to the earth 5 
months). Thus some longitudes were viewed more than once. The 
variation in the latitude of the subradar point during the entire 
period shown in Fig. 1 was within the range -3.6 to + 5.8° 

(assuming a pole position 6 - 66°, a = 270°). Since, however, 

50% of the echo power is returned in the first 300 y sec from 
the leading edge (Evans et al. 1965) corresponding to a region 
of ± 7° about the subradar point, to a first approximation we may 
ignore the variation of latitude with time and simply plot the cross 
section versus longitude. This is done in Fig, 2 employing a 
coordinate system proposed by I. I, Shapiro and colleagues. In 
this system, the prime meridian is defined as being the central 
meridian at 0 . Oh U. T, January 1, 1961. The pole position 
adopted is that given above, and a rotation period of -243.2 days 
has been assumed. Since the region from which the radar return 
chiefly arises is of the order of 15° in longitude, it is 
evident that much of the fine structure in this figure is 
spurious . 

Anomalous scattering regions on the surface of Venus have 
been detected by a number of workers , and their location studied 
extensively by Carpenter (1966). We have attempted to see if 
there is any correlation between the major changes in cross 
section shown in Fig. 2 with the passage of such features through 
the subradar region. This is somewhat complicated by the difference 
in the coordinate systems employed by Shaprio and Carpenter. 

Very crudely Carpenter’s longitudes will be obtained by sub- 
tracting 240° from those shown in Fig. 2. However, this 
transformation is not exact since Carpenter (1966) employed a 
pole position of (6 = 68°, a = 255°) and a rotation period of 
250 days. 

One suggested identification is that the dip in cross section 
at 162° longitude is associated with the passage of the feature 
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B of Carpenter's list through the subradar region. A second dip 
occurs when the feature G would lie on the central meridan, but 
if Carpenter's value for the latitude of G is correct, then it 
seems that this feature would lie outside the central region 
that chiefly governs the echo strength. The peak at 320° 

(labeled H) in Fig. 2 occurs at the same longitude as an increase 
at 3.8 cm wavelength was observed (next section). 

In sum, Venus rotates with respect to the Earth such that the 
longitude of the subradar point changes at a rate of between 
3° per day (at superior conjunction) and ^1° per day (near 
inferior conjunction). Some 50% of the echo power is returned 
from a region at the subradar point extending some 15° in 
longitude. Thus, we might expect significant changes in the cross 
section to occur on a time scale of the order of 5 days or 
longer. The accuracy and the sampling afforded by the published 
data is too poor to permit the construction of a reliable cross 
section versus longitude map. However, there is nothing in the 
data to suggest that what real variation does exist is not 
simply a consequence of changes in the roughness and/or reflect- 
ivity of the subradar region. 

Week to week changes of the roughness have been reported by 

Evans et al. (1965) and increases in roughness would be expected 

to lower the observed cross section. Thus major dips are 

consistent with the view that a large rough area has moved to 

occupy the subradar region. The converse probably does not 

hold. The average scattering properties of Venus are sufficiently 

near specular that it would not be possible to make a major 

increase in the cross section by placing a completely smooth 

region at the subradar point. This follows because the back- 

scattering gain g of a sphere with a gently undulating surface 
2 . 

is 1 + a where a xs the rms surface slope. For most regions 
on Venus a ^ 0.1 and thus the gain is indistinguishable from that 
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of a perfect sphere. It follows that increases in cross section 
must result from increases in the reflection coefficient at the 
subradar point (or of the roughness of a large area surrounding 
this region) . 


CENTIMETER WAVE OBSERVATIONS 

Whereas the cross section reported for meter and decimeter 

2 

observations lies in the range 0.1-0. 2 nr , Karp et al. (1964) 

2 . 

reported a cross section at 3.6 cm of 0.01 TTr . This was subsequently 
confirmed by Evans et al. (1966) who were able to show that, at 
least in part, this low value is caused by attenuation in the 
atmosphere of Venus. Observations at this wavelength thus hold 
special interest in that they may provide information concerning 
the nature of the atmosphere and possible day-to-day changes. 

Figure 3 shows the cross section measurements reported by 
Evans et al. (1966). These values are also listed in Table 2. 

Near close approach ( i . e . , up to day 50) the accuracy of 
measurement was limited by uncertainty in the performance of the 
radar (e.g.,the antenna gain, transmitter power, etc), together 
with the amount of atmospheric attentuation introduced by the 
earth's atmosphere. Subsequently, the signal-to-noise ratio 
deteriorated, and the uncertainty in the determination of the 
echo power increased as reflected by the error bars in Fig. 3. 

In what follows we discuss these three principle sources of error 
in order to determine whether the two large cross section 
increases observed (Fig. 3) may be attributed to errors of 
measurement . 
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Figure 3. Cross section measurements reported by Evans, et al. (1966) at 



TABLE 2 


C. W, CROSS SECTION MEASUREMENTS OF VENUS AT 3.8 cm 

(Evans et al . 1966) 


Date 


No. 

Echo Power 
(dbw) 

Flight Time 
(sec) 

Cross Section 
2 

db below Trr 

Cross Section 
(% irr^) 

Jan . 

18 

18 

-181.41 

274.85 

-18.41' 

1.44 

Feb. 

1 

32 

-182.65 

275.89 

-19,59 

1.10 

Feb. 

9 

40 

-183.48 

300. 54 

-18.73 

1.34 

Feb. 

15 

46 

-185 . 78 

308.77 

-19.56 

1 . 11 

Mar . 

21 

80 

-194.50 

557,45 

-19.12 

1.2 2 

Mar . 

30 

89 

-194.32 

626.59 

-16. 81 

2.08 

Apr . 

15 

105 

-197 .22 

752.89 

-15.34 

2.92 

May 

5 

125 

-202.70 

909.02 

-18.75 

1.33 

May 

17 

137 

-202.33 

1000.36 

-16.40 

2.29 

May 

28 

148 

-203.66 

1082.40 

-14, 84 

3.28 

June 

15 

166 

-207.09 

1208.53 

-18.07 

1.5 6 

June 

27 

178 

-208.72 

1285 . 56 

-18 . 88 

1.29 


Radar Calibration Errors 

Evans et al. (1966) estimated that the overall accuracy of 
the calibration of the radar is probably not better than ± 2 db. 
Thus , the absolute values reported may be in error by as much as 
± 50%. However, the operation of the equipment was monitored 
sufficiently well that variations in performance larger than 
± 0.5 db would not have escaped notice. Thus the relative 
accuracy of the measurements should not be worse than ± 10%. It 
follows that instrumental effects can be excluded as the prime 
cause of the variability shown in Fig. 3. 
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Signal-to-Noise Ratio Errors 

The spectrum of the signals (in the presence of noise) was 
determined by digital spectrum analysis of the output of a 500 
c/s wide channel (i.e., wider than the signal bandwidth) . The 
output of the same channel was also analyzed when the signal was 
absent, and this "noise only" spectrum subtracted from the 
"signal-plus-noise." To remove gain differences between the two 
runs , the two spectra were normalized before subtraction to have 
the same mean value in two frequency intervals located on either 
side of the band occupied by the echo. This operation introduces 
the largest uncertainty in the determination of the echo 
intensity. If the combined bandwidth of the two windows is b c/s, 
and the duration of a run is t seconds, there will be a baseline 
uncertainty after subtraction of AT where 

AT = T n /27bt °K (1) 

where T n is the system noise temperature (°K). If n runs are 
averaged, this becomes 

AT b T /2/nbt °K (2) 

n 

If the mean echo temperature over the frequency band occupied by 
the signal is T °K, then the uncertainty in the echo power is 
simply AT/T*. Table 3 lists values of T, T , b, nt, AT, and AT/T. 
It can be seen that the accuracy declined with time, although 
improvements in T^ and b were made which served to reduce AT. On 
June 15 and 27 the comparison bands were widened to the point 
where they occupied portions of the echo spectrum. Thus, the echo 
contributions from the limb region were not included in the 
computation of T, and the cross section values for these two days 
reported in Table 2 should be increased by about 10% to allow for 
this effect. 

*The uncertainty introduced by the stochastic nature of the signal 
itself has here been neglected. 


C-ll 



TABLE 3 


SIGNAL-TO-NOISE UNCERTAINTY 
IN THE 3.8 cm CROSS SECTION MEASUREMENTS 

(Evans et al. 1966) 


Date 

T 

(°K) 

T 

n 

(°K) 

b 

(c/s) 

nt 

( sec) 

/-N 

o 

* > 

w H 

AT/T 

(%) 

Jan. 18 

400.0 

156 

107.4 

3510 

0.36 

0.09 

Feb. 1 

325.6 

173 

107.4 

2700 

0.47 

0.14 

Feb. 9 

213.5 

185 

179.7 

42 0 0 

0.30 

0.14 

Feb. 15 

105.2 

237 

179.7 

2400 

0.51 

0.48 

Mar. 21 

8.89 

209.7 

93.7 

3240 

0.54 

6.0 

Mar. 30 

9.28 

119.2 

93.7 

5 6 7 0 

0.16 

2.1 

Apr. 15 

4.51 

123.1 

93.7 

5250 

0.25 

5.5 

May 5 

1.28 

121.6 

93.7 

5400 

0.24 

18.9 

May 17 

1.39 

127.4 

93.7 

8000 

0.21 

14.9 

May 28 

1.02 

197.5 

93.7 

7560 

0.33 

32.4 

June 15 

1.2 3 

126.5 

250 

7248 

0.14 

10.8 

June 27 

0.83 

114.0 

250 

6425 

0.09 

11.0 

Table 

3 shows 

that the 

uncertainty 

in the 

signal- 

-to-noise 


ratio cannot have contributed to the day-to-day variability to a 
greater extent than variations in system performance. These two 
sources combined would be expected to introduce an rms fluctuation 
of the order of 45% at most. This is far less than the fluctuation 
observed (Fig. 3). 
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CORRECTIONS FOR TERRESTRIAL ABSORPTION AND POINTING ERROR 


Two corrections were applied to the values of echo power 
(Table 2) in order to obtain the values for the cross section 
given. The first of these took account of the fact that the true 
and apparent positions of Venus differed. That is, the direction 
to which the antenna should be pointed to place the greatest 
amount of transmitter power on the surface is not precisely the 
same as that for receiving the largest echo signal owing to the 
motion of Venus. In the early measurements the antenna was 
maintained pointing in the apparent position, and although the 
apparent motion of Venus was greatest at this time, the flight 
time was sufficiently short that the pointing error so introduced 
was negligible. Later, by pointing at the apparent position the 
effective gain in the "transmit" interval was reduced. An 
estimate of this reduction was made from the antenna beam pattern 
and the angular separation between the true and apparent positions. 
The largest correction that was necessary to apply due to this 
effect was 0.5 db (Table 4) and the accuracy of this correction 
is judged to be better than £ 0.2 db. 

Atmospheric attenuation occurs in the earth's atmosphere even 
in fine weather. The one-way attenuation with the beam in the 
zenith is at most 0.1 db and at the zenith distance for which 
most of these measurements were made would be £ 0.2 db (Evans 
et al. 1966). Thus, the two-way loss would be £ 0.4 db. This 
correction was not applied to the results listed in Table 2. 

On some days observations were carried out during rain or heavy 
snow and the increased terrestrial atmospheric absorption present 
was reflected in higher system noise temperatures T . Corrections 
for the additional loss on these days were obtained as follows. 

The mean of the system temperature values encountered on fine days 
T n was obtained and the system temperature observed during bad 
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weather then used to compute the loss L from 

T sky (1 - L) + V + T n = T n t3) 

where T^y is the radio temperature of the sky background at this 
wavelength (taken to be 5°K) and is room temperature (taken to 

be 270°K). This is not a very satisfactory procedure and can be 
criticized on a number of counts. The proper way to proceed would 
be to determine the atmospheric extinction directly by observing 
the apparent cross section as a function of zenith distance. 

TABLE 4 

CORRECTIONS APPLIED TO THE 3.8 cm CROSS SECTION MEASUREMENTS 
FOR ATMOSPHERIC ATTENUATION AND POINTING ERROR 

(Evans et al. 1966) 


Date 


Pointing 

Correction 

Atmospheric 

Correction 

Total 

Correction 

Jan . 

18 

0 

db 

0 

db 

0 

db 

Feb. 

1 

0 

db 

0 

db 

0 

db 

Feb. 

9 

0 

db 

0.2 

db 

0.2 

db 

Feb. 

15 

0 

db 

1.2 

db 

1.2 

db 

Mar . 

21 

0 .1 

db 

0 

db 

0.1 

db 

Mar . 

30 

0.2 

db 

0 

db 

0.2 

db 

Apr . 

15 

0.2 

db 

0 

db 

0.2 

db 

May 

5 

0 

db 

0 

db 

0 

db 

May 

17 

0.45 

db 

0 

db 

0.45 

db 

May 

28 

0.5 

db 

1.3 

db 

1.8 

db 

June 

15 

0 

db 

0 

db 

0 

db 

June 

27 

0 

db 

0 

db 

0 

db 
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An examination of the data has been undertaken in order to 
attempt this. Unfortunately, the signal- to-noise ratio of the runs 
taken beyond March 30 is too poor to permit one to detect small 
changes in cross section from run to run. Prior to this date the 
elevation of Venus was low (^30°) and data taking generally 
stopped when the elevation fell to 20°, Thus, the range of 
zenith distances over which observations were made is small. A 
further difficulty is that the system temperature T n varied with 
zenith distance, but was not usually measured at frequent enough 
intervals to permit the cross section to be determined accurately 
for individual runs. (This follows because the value of T 

n 

assumed for the baseline determines the overall calibration of 

the signal spectrum.) Figure 4 shows an attempt to obtain a 

relation between the echo temperature and the zenith distance for 

February 9. Two sets of points are shown. The open circles 

represent values obtained by assuming a constant value of T n 

(taken to be the mean of all the measured values)- The closed 

circles (thought to be more accurate) represent values obtained 

by taking estimates of T for each run (obtained by reasonable 

n 

extrapolation through the measured values). A similar attempt with 
the data of March 30 is presented in Fig. 5. The straight lines 
shown in Figs. 4 and 5 have been fitted by eye. That shown in 
Fig. 4 implies a two-way zenithal absorption of 0.44 db and 
Fig. 5 some 0.56 db. The uncertainty in the value derived for 
March 30 is considered far larger than that for February 9 in 
view of the wider scatter of the points (Fig. 5), 

R. Allen (1967) has examined the variation of the intensity 
of a number of radio sources observed with the Haystack antenna 
as a function of zenith distance x* Each measurement was corrected 
for the atmospheric extinction according to a theoretically 
obtained dependence. 

F = F q exp (-.025 sec %) (4) 
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Figure 4 Variation of the equivalent echo temperature of Venus at 3.8 cm with z 
angle x on Feb. 9, 1966. Open circles represent values obtained assuming a cons 
receiver temperature. Closed circles are values based upon estimates of how the 
receiver temperature varied with zenith distance and should be more accurate. 
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Figure 5 Variation of the equivalent echo temperature of 
Venus at 3.8 cm with zenith angle on March 30, 1966. 





where F is the observed flux and F would be the flux in the 

o 

absence of atmospheric attenuation. According to this law the one- 
way zenithal abosrption is 0.055 db. When the corrected values 
of flux were plotted against x a residual variation was found 
which was attributed to variation in the antenna gain. Figure 6 
shows the dependence in the observed gain G as a function of the 
elevation angle between 20° and 40°. Allen’s results imply that 
over this elevation range 

G - G q exp (-.0256 sec x) (5) 

The similar dependence of this expression to that for 
atmospheric extinction (Eq. 4) suggests that the latter effect 
may have been underestimated. 

The antenna gain employed in the radar equation in order to 
reduce the observed echo power to a radar cross section was 
established from observations at high elevations. Thus these two 
expressions should be combined, and allowing for the two-way 
nature of both effects we would expect the radar cross section to 
vary as 

a = £exp (-0.101 sec x + 0.051)j (6) 

This is close to the law determined empirically (Figs. 4 and 5). 

Since most of the radar observations of Venus were conducted at a 
zenith distance of ^ 60° (sec x = 2) we believe that the values 
reported in Table 2 are systematically low by approximately 
0.8 db. 

Adopting an average value of 0.11 db for the one-way attenuation 
we can estimate that the atmosphere contributes some 6.7° to the 
system temperature T n even under fine conditions. This is suf- 
ficiently small that the expression employed to derive the 
correction for atmospheric extinction (Eq. 3) in wet weather is 
probably not seriously in error. In retrospect, however, it seems 
unwise to apply large corrections to the result in view of the 
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ELEVATION (deg) 


Figure 6 The variation of the gain of the Haystack antenna with 
elevation angle according to radio star observations (Allen, 1967). 
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cumulative effect of errors in T . If by chance an erroneously 

high value of T is measured, then a cross section will be derived 

n 

which is correspondingly high. This will then be scaled by a factor 
to allow for atmospheric attenuation which is also an overestimate. 
Accordingly, it seems best to discard the values reported for 
February 15 and May 28, and Table 5 presents a revised list of cross 
sections in which this has been done. In this table a correction 
of 0.8 db has uniformly been applied to the- results to allow for 
the combined effects of atmospheric extinction and gain degradation. 

TABLE 5 

REVISED VALUES FOR THE CROSS SECTION AT 3.8 cm 


Date 


Uncorrected 
Value 2 

(db below fir ; 

Atmospheric 
Correction 
i (db) 

Corrected 

Value 

(db below nr 

2 > 

(% nr 2 ) 

Jan. 

18 

-18.41 

-0.8 


-17.61 


1.74 

Feb. 

1 

-19.59 

-0.8 


-18.79 


1. 32 

Feb. 

9 

-18.73 

-0.8 


-17.93 


1.61 

Mar . 

21 

-18 . 32 (1) 

-0.8 


-17.52 


1.77 

Mar . 

30 

-16.81 

-0.8 


-16.01 


2.51 

Apr . 

15 

-15.34 

-0.8 


-14.54 


3.50 

May 

5 

-18.75 

-0.8 


-17.95 


1.60 

May 

17 

-16.40 

-0.8 


-15.60 


2.75 

June 

15 

-17. 57 (2) 

-0.8 


-16.77 


2.10 

June 

27 

-18. 38 (2) 

-0.8 


-17.58 


1.75 


(1) 

Revised value; 

allows for 0.8 

db 

estimated 

additional 



waveguide loss 

on this day. 






(2) 

Revised value; 

allows for 0 . 5 

db 

reduction 

in 

cross 



section due to 

narrow bandwidth. 
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Longitude Variation 


Table 5 contains three large values (March 30, April 15, and 
May 17). The mean of the remainder is 0.017 Trr . Thus the 
largest increase is of 0 the order of 100%. Similar increases 
have been observed at decimeter wavelengths (Figs, 1 and 2). 

Fig. 7 compares some of the values listed in Table 5 with the 
cross section observed at 23 and 70 cm (Fig. 2). The data are 
too sparse to draw any firm conclusions, but we do notice that 
the first increase in the 3.8 cm cross section occurred roughly 
at the longitude for which the decimeter cross section appears to 
be highest. In conclusion, we believe that cross section 
variations of a factor of two do occur' at 3.8 cm wavelength, but 
it is impossible to decide from the foregoing whether these result 
from changes in the absorption by the atmosphere of Venus or from 
variations in the surface reflectivity. 

A MAJOR SURFACE FEATURE ON VENUS 

The spectra obtained for Venus near inferior conjunction by 
Evans et al. (1966) at 3.8 cm were carefully examined in order to 
obtain the mean scattering function. Subsequently the spectra 
were not examined in any great detail though Evans et al. (1966) 
observed a tendency for them to become assymetrical prior to an 
increase in the cross section. 

Fig. 8 shows the average echo spectrum observed on March 21 
together with the mean curve expected based upon earlier obser- 
vations. Difficulties were encountered with the doppler tuning 
on this day, and the echo power to the right of the expected 
curve was thought to result from instrumental smearing. Figure 
9 shows a spectrum obtained for the same day by combining only 
runs during which the doppler compensation was known to have been 
applied correctly. A large feature can be seen on the approaching 
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VENUS CROSS SECTION 







ECHO 



Figure 9 Spectrum of Venus signals March 21, 1966, with runs 
in which doppler tracking was believed poor removed. 
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limb and one on the receding limb. The latter can be identified 
with G of Carpenter's list, whife the former has here been labeled 
H and I at the points of greatest intensity. On March 30 (Fig. 10) 
the approaching feature had moved to within 10° of the subradar 
point and extended westwards to at least 20° in longitude and 
possibly further. By April 15 (Fig. 11) the identifiable peaks H 
and I lay on the receding limb. Based upon the assumption that 
the regions responsible for these reflections lie close to the 
equator, we find that they move in the expected manner. The 
longitudes obtained for the two peaks are H = 332 ±1°; I = 322 
± 1 °. 

It is possible to compute what contribution these features 

make to the total cross section assuming that the area beneath 

2 

the broken curve corresponds to the mean value of . 017irr . The 
increases are, respectively, 12, 36 and 35% for March 21, 30 and 
April 15. The cross sections measured (Table 5) correspond to 
increases of 5, 60, and 100% above the mean. In the case of 
April 15 the proper position of the expected curve is by no 
means obvious, and it is possible that it should be lower than 
shown, thereby reducing the above discrepancy. Such might be the 
case, if (as seems likely) the feature extends a further 20° 
westward beyond I. In this event the feature would occupy a 
40° region in longitude and be considerably longer even than the 
complex region g discovered by Goldstein (1965). 

In the case of March 30 the difference in the cross section 
derived directly and by comparison with the mean curve is 0.7 db 
which is certainly within the experimental accuracy expected. 

In sum, it seems that the increase in cross section observed 
early in April is associated with the appearance of a large 
feature on the disc of Venus. When sufficiently removed from 
the region of the subradar point the increase in cross section 
is probably attributable to the increased power contributed by 
this feature. It follows that we cannot hope for perfect 
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correlation between the results of the decimeter and centimeter 
measurements compared in Fig. 7. In the pulse-type decimeter 
measurements, the amplitude of the echo attributable to the near 
surface was taken as a measure of the overall cross section, it 
being assumed that the scattering law does not change with time. 

It is interesting that the highest cross section thus far 
reported at decimeter wavelengths occurred when the peak H lay on 
the meridian (Fig. 7). Thus it must be supposed that the region H 
has a higher reflectivity than average. Yet we are also required 
to make region H rougher than average in order to account for why — 
relative to its environs — it gets brighter the further it lies from 
the center of the disc (Figs. 9-11). The increase in cross section 
observable at 3.8 cm wavelength may result from reflectivity and 
roughness variations and/or height variations of the surface 
resulting from changes in the atmospheric absorption. Precise 
cross section measurements at two wavelengths together with 
accurate range measurements should help to resolve this question. 

SUMMARY 

Variations of the cross section of Venus have been reported by 
a number of observers (Table 1). Because of its long rotation 
period with respect to the earth (^ 5 months) it has been difficult 
to construct a reflectivity vs longitude plot for Venus corresponding 
to those now available for Mars. Difficulties arise in maintaining 
the performance of the radar constant (or at least known) and from 
the large variation in the absolute signal strength due to range 
changes. Figure 2 represents a first crude attempt at such a 
plot and includes data obtained at 23 and 70 cm. Much of the fine 
structure in Fig. 2 must be spurious since the region from which 
50% of the echo power is returned occupies some 15° in longitude, 
and thus there should be good correlation over longitude intervals 
of half this amount. 
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Variations in the cross section are to be expected depending 
upon the roughness of the terrain occupying the subradar region. 
Direct observations of roughness changes are available from 
scattering law measurements (Evans et al. 1965) and from the 
detection of anomalous scattering regions on the surface 
(Carpenter, 1966), some of which appear to move through the 
subradar point. Increasing the roughness at the subradar point 
can readily account for a reduction in the cross section by a 
factor of 2 . An increase in the cross section must largely be 
caused by increased reflectivity. 

Measurements at 3.8 cm by Evans et al. (1966) have been re- 
examined in this paper. The possible sources of error in these 
measurements were discussed and it was concluded that attenuation 
in the earth's atmosphere and degradation of the antenna gain with 
zenith distance were the most serious of these. As a result, two 
of the measurements were discarded and the remainder were revised 

upward (Table 5). The revised values exhibit increases above the 
2 

mean (0.017 irr ) of the order of 100%. By examining the spectra 
of the signals obtained during one such instance, it was found 
that the increase appeared to be associated with the passage across 
the disk of a large anomalously scattering region, occupying 
approximately the longitude interval 320° - 350° + 10°. We are 
led to conclude, therefore, that fluctuations in the cross section 
at this wavelength are probably related to variations in the nature 
of the terrain visible to the radar (e.g. , reflectivity, roughness, 
and perhaps height) , and are not the result of weather-like 
phenomena in the atmosphere. 
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RADAR OBSERVATIONS OF VENUS AND MARS 


ABSTRACT 

Planetary radar observations commenced in 1960 with the suc- 
cessful detection of echdes from Venus by groups in the U. S., 

U. K. , and the U. S. S. R. Since then Venus has been extensively 
studied by radar, and similar though less extensive measurements 
have been made on Mercury and Mars . 

Possibly of greatest utility to the space program are the 
revisions obtained by radar to the values for the astronomical 
unit, the earth-moon mass ratio, the radii of the planets and the 
elements of their orbits. These revisions make it possible to 
direct a deep-space probe to within the immediate vicinity of the 
planet where earlier optical data would have been grossly in error. 

In addition to the improvements to planetary positions radar 
has provided new values for the rotation periods of Mercury and 
Venus, and estimates of their surface roughness and composition. 

In the case of Venus some information concerning the absorption 
of radio waves by the atmosphere has also been obtained. For Mars 
the reflectivity and surface elevation vs Martian longitude 
have been determined for a region in the Martian tropics. These 
quantities show interesting correlation with visual (i.e. dark) 
features. 
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RADAR OBSERVATIONS OF VENUS AND MARS 
J. V. Evans* 

Lincoln Laboratory^ 

Massachusetts Institute of Technology 
July 1967 

INTRODUCTION 

It is the purpose of this paper to review the status of 
recent radar observations of the planets Venus and Mars and to 
point out some of the logical extensions of these radar inves- 
tigations. Further experimental data and its interpretation for 
the planet Venus is contained in a companion paper (See Appendix C 
Evans, J. V., "Variations in the radar cross section of Venus" 

July 1967). The support of the staff of the Haystack Microwave 
Facility in the work described here is gratefully acknowledged. 

We are also indebted to Dr. Irwin Shapiro of Lincoln Laboratory 
whose staff kindly supplied us with the working ephemerides 
necessary for the taking of radar data. The work of preparing 
this review was performed during the course of the "TYCHO" Study 
Group meeting, Dartmouth College, Hanover, New Hampshire, June - 
July, 1967. 

RADAR CONTRIBUTIONS TO PLANETARY SIZES AND MOTIONS 

Size and Distance 

Precise radar distance and velocity measurements have been 
made by a number of observers (e.g. Evans et al, 1965, 1966a 
Pettengill et al. 1962, 1967). These have been reduced in con- 
junction with appropriate optical observations of the planets 
(by the U. S. Naval Observatory) to yield revisions for the 
orbital elements of Earth, Venus and Mercury, the Earth-Moon mass 
ratio, the radii of Venus and Mercury and the astronomical unit 

* Presently Visiting Professor, University of Illinois, Urbana, 111. 

# Operated with support from the U. S. Air Force. 



(Ash, Shapiro, and Smith 1967). 

The methods of making precise delay and velocity measurements 
have been well documented by the above observers, as has the pro- 
cedure by which these data have been reduced. Thus, we will 
merely summarize the conclusions reached in the form of Tables. 

Table 1 presents the most recent value for the astronomical unit 
(presented in light seconds) together with values for the radii 
of the planets Mercury and Venus. 

TABLE 1 (After Ash et al. 1967) 

THE ASTRONOMICAL UNIT AND PLANETARY RADII 

General Formal Standard 

Newtonian Relativity Error 

A. U. (light-sec) 499.004785 499.004786 5 x 10~ 6 

Mercury Radius 2440.0 2434.0 2.2 

(km) 

Venus Radius 6055.5 6055.8 1.2 

(km) 

Two estimates are presented for two hypotheses concerning the nature 
of the universe. In the first it is assumed that the motion of 
the planets proceeds according to the laws of Newtonian gravitation, 
and that light propagates always as in flat space, whereas in the 
second the consequences of general relativity on both the motion 
of the planets and the interaction of the electromagnetic probing 
signal with the gravitational field of the sun are assumed to hold. 
Table 2 presents the inverse masses of the inner planets and the 
earth-moon mass ratio also based upon the combined radar and optical 
data as processed by Ash et al (1967). In both tables the formal 
standard errors quoted represent the internal consistency of the 
results, and undetected systematic errors may exist which are 
significantly larger. 
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TABLE 2 (After Ash et al. 1967) 


INVERSE MASSES 

OF INNER PLANETS 

AND EARTH-MOON 

MASS RATIO 


Newtonian 

General 

Relativity 

Formal Standard 
Error 

Mercury 

6029000 

6021000 

53000 

Venus 

408450 

408250 

120 

Earth and Moon 

328950 

328900 

60 

Mars 

3106700 

3111200 

9 000 

Earth-Moon mass 
ratio 

81.3024 

81. 3030 

0.005 


Orbital Elements 

Although 15 years of optical data were employed in this study 
(Ash et al, 1967) and only a few years of radar data, the latter 
serve to improve the initial mean anomaly, and the orbital 
eccentricity, and to redetermine the semimajor axis for each planet 
with comparable accuracy. These resuls are summarized in Table 3, 
where only the "relativity fit" is given. Since only a few pre- 
cise distance determinations of Mars have been made thus far, no 
results have been published for radar/optical determinations of 
the orbital elements of that planet. 

Rotation Rate and Pole Position 

The rotation rate of Mars is well known from optical 
observations and radar is unlikely to make any significant improve- 
ment. In the case of Venus, however, radar provides the only 
reliable Earth-based technique for determining the sense and rate 
of rotation. For this purpose several different radar techniques 
have been employed. Goldstein (1964) and Carpenter (1964) 
attempted to recognize the echo power in the frequency spectrum 
associated with the planet’s limbs and thus obtain directly the 
maximum doppler broadening. This so-called "base-bandwidth" 
technique led to a value of 250 40 days retrograde, which was 

later refined to 250 +_ 9 days (Goldstein 1965a). A second method 
which takes advantage of radar "features" which are seen in the 
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OSCULATING ELLIPTIC ORBITAL ELEMENTS AT JED 2439340. 
REFERRED TO MEAN EQUINOX AND EQUATOR OF 1950.0 
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doppler broadened power spectra and which correspond to regions on 
the surface that reflect anomalously. Because Venus presents 
approximately the same face towards Earth at successive inferior 
conjunctions, these tend to reappear from conjunction to conjunction. 
Employing this "feature method" Carpenter (1966) concluded that the 
period lay between 244 and 254 days. Subsequently, Evans et al. 
(1966b) and Carpenter employed the method over two successive 
conjunctions to obtain a period of 244 days, to within an error 
of only about 1%. The most recent determination reported in the 
literature (Dyce et al, 1967a) is based upon combined delay- 
doppler mapping and yields a value of 244.3 +_ 2 days (retrograde). 

The pole position is at declination -■6.6,4° +_ 1° and right ascension 
90.9° +_ 1° (1960 epoch). The inclination of the axis is about 87° 
with respect to the orbital plane of Venus and approximately 90° 
with respect to the plane of the Earth’s orbit (i.e. , the 
ecliptic) . The coincidence of these results (i.e. , the axis 
inclination and rotational period) with the theoretical orientation 
and period (of 243.16 days) which would result were the Earth's 
orbit to control the rotation of Venus through an interaction 
with a permanent dipole moment of that planet is remarkable. In 
the event this locking proves to exist, Venus would, on average, 
just complete 4 axial rotations between inferior conjunctions as 
seen by a terrestrial observer. The length of a Venus "day," 
i.e., the rotation period with respect to the Sun, corresponds to 
about 177 Earth days. 

The rotation period of Mercury was long thought to be close 
to 88 days, i.e., equal to its orbital period. Delay-doppler 
measurements (Pettengill and Dyce, 1965), however, showed that the 
period is more nearly 59 (+ 3) days, and that Mercury, therefore, 
executes 3/2 axial rotations in each orbital revolution (Colombo 
and Shapiro, 1965). This revision is now supported by recent 
optical evidence (see Dyce et al. , 1967a). All that can be said 
concerning the tilt of the pole is that it appears to be normal 
to the plane of Mercury’s orbit to within 62°. If the planet is 


D-5 



truly locked, as seems highly probable, the solar day on Mercury 
would be 176 Earth days (i.e. , just twice the orbital period, 
since the 3/2 resonance requires that the planet present alternate 
faces to the sun at successive perihelion passages). 

RADAR CONTRIBUTIONS TO PLANETARY SURFACES 

Venus 

Surface material . - The radar cross section of Venus has been 
obtained at a number of wavelengths and is plotted in Fig. 1 
(Evans et al. 1966b). The dotted curve has merely been put in by 
eye. For comparison the moon exhibits no clear wavelength 
dependence but has a cross section a = 7% over the range 
1 cm < X < 10 m. It is tempting to conclude that the marked 
dependence shown in Fig. 1 results from atmospheric attenuation 
(next section) though it is not impossible to construct a surface 
having reflection properties which would give rise to the reflec- 
tion law indicated. Such a surface would, however, require care- 
ful control of the dielectric constant (both real and imaginery 
parts) as a function of depth, and would require a remarkable 
degree of horizontal uniformity, such that it seems it can be 
discounted as being exceedingly artificial. We cannot discount, 
however, the possibility that the trend shown in Fig. 1 results 
from a combination of atmospheric attenuation and reflection co- 
efficient changes occurring with wavelength. 

The mean of the long wave (X > 23 cm) measurements is a = 16%, 
which corresponds to a dielectric constant of e ^ 5.5. This value 
would be consistent with dry compact rocklike material and definitely 
excludes water as the surface covering. It is interesting to note 
that the reflection coefficients of the Moon, Mercury and Mars on 
average lie between 6% and 8%. It seems unlikely (though of course 
not impossible) that the surface material on Venus is entirely 
unlike that on the other three bodies. A more plausible argument 
is that Venus is distinguished by having more compact surface 
material than the other three. That is, it may be supposed that 
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the surfaces of the Moon, Mercury, and Mars have been eroded (we 
have evidence of this in the case of the Moon and Mars) such that 
the surface rocks have been ground to particulate matter with a 
porosity of 50% or more for some considerable depth (centimeters). 

A collorary is that on Venus such processes are not at work or that 
the surface rocks melted and resolidified in the not too distant 
past. 

Surface roughness .- The angular scattering law for Venus has 
been determined in different ways by a number of observers. As an 
example, Fig. 2 shows that the mean echo power vs delay obtained 
by Evans et al. (1966a). The measurements show a fall off in echo 
power with delay which is more rapid than that from the moon. 

Since the delay is related to the cosine of the angle of incidence 
<j> , it is possible to replot Fig. 2 as echo power vs <p and regard 
such a plot as a spectrum of surface slopes (at least for c|> < 30°). 
On this basis the mean slope obtained for Venus when averaging 
over-all azimuthal planes is 8.2°. The corresponding value for the 
moon at X = 23 cm is 10.2°. These slopes probably correspond to 
intervals over the surface of the order of 10 X, or say 2 meters. 
When sampled over smaller intervals, the mean surface slope for 
the moon is found to increase and conversely when the wavelength 
is increased the mean slope decreases. No such clear-cut effect 
exists for Venus. Fig. 3 presents echo power vs delay curves 
similar to that of Fig. 2 obtained at 70 cm (Pettengill, 1965), 

23 cm (Evans et al. 1965) and 12.5 cm (Muhleman, 1965). The 
observations at 70 cm and 23 cm were made directly using short 
pulses, whereas the 12.5 cm curve was obtained from the mean echo 
power frequency spectrum by an appropriate transformation (Evans 
et al. 1965). This difference, together with variations in the 
pulse lengths employed (in the case of the pulse radars) caused 
the resolution achieved in the region <_ 1.0 m sec delay to differ 
between the three experiments. This probably accounts for most 
of the divergence between the curves near the origin (Fig. 3). 
Another factor that enters is that the scattering law actually 
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Figure 3 . Comparison of the echo power vs delay for Venus at three 
wavelengths. The divergence of the curves for delays of less than 1 
msec is believed in part to be due to differences in the range 
resolution achieved at the three wavelengths and in part to differences 
in the nature of the terrain near the subradar point when the three 
sets of measurements were made (Evans, et al. 1966b). 




changes with time depending upon the nature of the terrain at the 
subradar point, and none of the curves shown in Fig. 3 represents 
a long-term average. Thus, it is concluded that no systematic 
variation of scattering has yet been established for Venus, and hence 
the spectrum of surface roughness (or the number or different types 
of surface feature) must be less than for the Moon. 

At large angles (<t> > 6 0°) the echo power has been attributed 
solely to surface structure having horizontal and vertical dimen- 
sions comparable with the wavelength. As such this structure can 
scatter into wide angles and provides the predominant scattering 
mechanism at grazing angles, because then so few large surface 
elements are tilted normal to the ray. In the case of Venus at 
23 cm wavelength about 11% of the echo power is associated with 
these scatterers. For the moon at the same wavelength the corres- 
ponding figure is twice as much. Thus not merely are the surface 
elements of Venus more gently undulating than those on the Moon, 
they are covered to a lesser extent by boulders , small craters or 
comparable structure. 

Surface features. - Regions on the moon — such as the rayed 
craters — are known to be anomalously scattering in that they are 
considerably brighter than their environs. Similar regions, though 
of unknown character, are visible in the radar returns from Venus. 
Figure 4 shows three such regions visible in an echo power vs delay 
plot such as was presented in Fig. 2. The greatest amount of 
attention to these features has been provided by the Jet Propulsion 
Laboratory (J.P.L.) group (e.g. Carpenter 1966) who have observed 
them as departures from the mean echo power vs frequency spectra 
(e.g. Fig. 5). In the early J.P.L. observations (Goldstein 1965a) 
only two features (named a and 6) were resolved. Carpenter (1966) 
subsequently showed that the region 6 was complex and contained at 
least 3 revolable features (seen on the lefthand side of Fig. 5). 

In all. Carpenter was able to identify some seven scattering centers 
and by plotting their doppler shifts vs time (Fig. 6) was able to 
locate them on the planetary surface (albeit with ambiguity in some 
cases). Table 4 presents the positions deduced for these features in 
a coordinate system defined by Goldstein (1965a) . 
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Figure 6 . The doppler vs date histories of the features recog- 
nized by Carpenter (1966) and listed in Table 4. Reading from 
left to right at about January 18 these are A through G. Ambiguity 
is introduced by the blanding of the features in early January which 
prevents precise locations being assigned to features B through D 
t Carpenter 1966). 
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TABLE 4 (After Carpenter 1966) 

LOCATION OF ANOMALOUS SCATTERERS ON THE SURFACE OF VENUS 


Feature Designation 


Longitude 


Latitude 


Goldstein 

(1965) 


Carpenter 

(1966) 



A 

-108.6 

+ 

.44 

+ 6.3 

+ 

2.3 

B 

- 78.1 

+ 

2.6 

- 3,3 

+ 

8.7 

B1 

- 75,8 

+ 

.62 

+ 11.9 

+ 

4.4 

C 

- 68.9 

+ 

1.3 

- 6.8 

+ 

5.8 

Cl 

- 76.8 

+ 

.63 

+ 16.7 

+ 

2,1 

C2 

- 70.0 

+ 

.65 

- 2,5 

+ 

3. 7 

D 

- 61.0 

+ 

1.8 

- 1.7 

+ 

7.1 

D2 

- 70.0 

+ 

.65 

+ 22.7 

+ 

1.7 

D3 

- 60.5 

+ 

,67 

- 2.5 

+ 

4.1 

E 

- 49.3 

+ 

.92 

+ 6.4 

+ 

8. 5 

F 

00.0 

+ 

.69 

-26.7 

+ 

1.8 

G 

+ 15.7 

+ 

,83 

-12.4 


2.1 


In this coordinate system the prime meridian (zero longitude) 
was defined as the central meridian as seen on 23 July 1964, 
when it then passed through a. The features seen in Fig. 4 are 
believed to be (from left to. right) F, G , and A. Those in 
Fig. 5 are B, C, D, and F (again reading from left to right). 

The feature to the right of F in Fig. 5 that is not marked with an 
arrow is thought to be G. It is interesting that the features are more 
readily observable in the depolarized echo component. This 
indicates that (like Tycho in the case of the moon) they cannot 
be explained as merely regions of anomalously high reflecting 
material but must also be rougher. That they seem more prominent 
in the 3.8 cm wavelength spectrum than 12.5 cm (Fig. 5) would 
be consistent with this , but one could also explain this as a 
consequence of a height difference between these regions and 
their environs. If the anomalous regions were mountains then 
the lower atmospheric absorption to the mountain tops would 
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render them brighter than their environs by a larger amount at 
3.8 cm than at 12.5 cm where little absorption may be presumed to 
take place. 

The similarity of the spectra shown in Fig. 5 is one argument 
to discount the possibility that the 3.8 cm echo arises from 
scatterers in the atmosphere. Since features are seen in the 
spectra at both wavelengths and these move across the disc as the 
planet rotates we have strong evidence that both reflections are 
from a solid surface. 

In Fig. 5 the dotted lines represent the locus of the doppler 
shift (measured with respect to the center of the spectrum) of 
each feature as a function of date. Reading across at Jan. 18 
the lines correspond to the features A through G (Table 4). 

The features blend between January 10 and 20 making it impossible 
to determine the locations of B through D unambiguously and several 
alternate positions are listed in Table 4. Similar work at the 
Arecibo Ionospheric Observatory (A.I.O.) (Dyce et al. 1967a) 
suggests that the feature a (Feature F of Table 4) is probably 
less than 900 Km in extent E-W and in excess of 3800 Km N-S 
and confirms that region 3 is complex. 

In the companion paper the discovery of a major new feature 
is presented which is larger in longitudinal extent (^ 40°) 
than any hitherto observed. This feature lies approximately 
60° to the west of feature G (i.e. at about +76° longitude) and 
hence is not observable until about 50 days following inferior 
conjunction. As such it has not been reported previously. This 
new result depends upon a more careful examination of spectra 
obtained by Evans et al. (1966b) and 3.8 cm wavelength using the 
Massachusetts Institute of Technology Lincoln Laboratory Haystack 
microwave facility. In the case of this extended feature there 
is reason to believe that it is elevated with respect to the 
surrounding terrain. 
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Mars 


Surface material .- Mars has been observed by radar at a 
variety of wavelengths since 1963: 70 cm by Dyce et al . (1967b), 

40 cm by Kotelnikov et al. (1963), 23 cm by Evans et al . (1965), 

12.5 cm by Goldstein and Gillmore (1963) and by Goldstein (1965) 

and at 3.8 cm by G. H . Pettengill (private communication). Unfortunately, 

in some of these detections the signal to noise ratio was too 

poor to permit a reliable cross section determination and as yet 

it is not possible to construct a meaningful cross section 

vs wavelength curve (cf. Fig 1 for Venus) for Mars, As far 

as can be determined the average cross section is of the order 

of 7 - 81 at all wavelengths--a result akin to that for 

the moon. This cross section implies a dielectric constant of 

the order of 3 suggesting that the Martian soil is unconsolidated 

and fairly dry. 

Variations of the reflection coefficient are found at 
different Martian longitudes (Goldstein and Gillmore 1963, 

Goldstein 1965b and Dyce et al. 1967b). Fig. 7 shows a com- 
parison of the reflectivity observed by the Arecibo (70 cm) 
and J.P.L. (12.5 cm) groups respectively during the 1965 opposition. 

The J.P.L. data have been averaged over 10° intervals of longi- 
tude, whereas the A.I.O. data have been plotted individually and 
then the shaded region drawn to encompass all the points. 

Fig. 7 shows that there are variations in cross section 
of the order of 3:1. The lowest values (^3%) are lower than 
encountered for the Moon while the highest (^13%) are about 
twice that of the Moon and approach the high average value for 
Venus discussed above. The exact reason for these variations 
is not known, but the correlation of the peaks with the dark 
markings (Fig. 6) is conducive to speculation (e.g. Sagan et al. 

1967). The different optical albedo suggests changes in surface 
material and hence in the intrinsic reflection coefficient. 

However, changes in surface roughness are also believed 
responsible (see below). 


D-17 



MARE TRIVIUM SYRTIS 

ACIDALIUM CHARONTIS MAJOR 



( JO !N30U3d ) f>d 


D-18 


Figure 7 . The variation of the radar cross section of Mars at 70 cm (Arecibo 
Ionospheric Observatory, shaded region) and at 12.5 cm (Jet Propulsion Laboratory, 
dotted line) during the opposition of 1965 as a function of the central meridian 
of the visible disk. A strip map appears across the top showing the principle 

optical features of Mars near the sub-radar latitude of 20°N. (After Dyce , Pettengill 
ianchez, 1967) ■ 




Surface roughness .- The average angular scattering law (cf. 
Fig. 2 for Venus) has not yet been obtained for Mars. The smaller 
size and higher spin rate of Mars (compared with Venus) render it 
a much more difficult object to study. Thus there has been 
adequate signal-to-noise ratio in existing experiments to examine 
the power spectrum (cf. Fig. 5 for Venus) of the echoes only 
around the peak. It is found that the spectral width is narrowest 
when the echo intensity is high suggesting that, in part, the. 
increase in cross section arises as the result of an unusually 
smooth region occupying the subradar point. At these times the 
width of the spectrum indicates that the rms surface slope is of 
the order of only one or two degrees. This suggests reflection 
from a surface that has been leveled in some fashion as a salt- 
flat has on earth. 

The average spectrum (i.e. averaging over all longitudes) is 
somewhat broader, but still indicates that the surface of Mars 
contains fewer elements than the Moon with slopes >_ 5°. Beyond 
this qualitative statement little else can be said at present, and 
we need considerably greater sensitivity in order to measure the 
angular scattering law over a wider range of angles. 

Surface features .- Beginning in early April, 1967, short-pulse 
(60 microsecond) echoes were obtained from Mars using the Haystack 
microwave facility at 3.8 cm wavelength (G. H. Pettengill, private 
communication) . With this system direct measurement of the 
Martian topography at about 21°North latitude was possible as the 
planet rotated under the radar beam. The planetary rotation was 
sufficiently rapid that errors in the planetary ephemeris did not 
affect the determination of the relative heights of the distant 
surface. In a given night about 100° of longitude could be 
observed in this way. Observations spaced a week apart gave an 
overlap of about 40° and in every case the results from the over- 
lapping region for the two nights agreed well. 

Although the resolution offered by the use of a pulse length 
of 60 microseconds might appear at first glance to be only about 
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8 km, in most cases the bulk of the echo could be shown (from 
the narrowness of the early maximum) to occupy a very much smaller 
interval. The accuracy of the measurement was better than 1 km in 
almost all cases. Both the accuracy and resolution are obviously 
better where the signal was strong. 

Figure 8 shows the results obtained (G. H. Pettengill private 
communication). The heights plotted are relative, of course, and 
the zero height is arbitrary. Eventual reductions involving an 
improvement in the full orbit using the methods described by Ash 
et al. (1967) should yield an accurate radius for Mars at the zero 
reference level shown here. The relative intensity of the initial 
echo maximum, whose location in delay yields the height information, 
is also indicated on the diagram by the thickness of the line 
traced through the observational points. The locations of the 
more prominent dark features at the observed latitude (21°N) are 
given as well as some of the traditional names for others. 

Perhaps the most obvious conclusion that can be drawn 
immediately is that the dark regions do not correspond uniquely 
to elevated nor to depressed areas, at least at this latitude. 

This is in contrast to the conclusions reached by Sagan et al. 

(1967) on the basis of the J.P.L. cross section and spectral 
results alone. In their model the dark regions of Mars are 
highland areas. Such a model was originally proposed by Wells 
(1965) on the assumption that the lighter areas of Mars result 
from a layer of fine dust and that these particles preferentially 
settle in the lowlands. Reasoning from the slight displacements 
between the longitudes of the radar maxima and nearby visually 
dark regions, Sagan et al. (1967) conclude that the dark areas 
correspond to regions of the Martian surface which are elevated 
by some 5 to 15 km above the lighter areas. Such is not the case 
apparently. 

It is also evident from Fig. 8 that the signal maxima occur 
not precisely but only approximately when the dark areas 
lie under the subradar point. The steepest slopes that may be 
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derived from the data in Fig. 8 are in the vicinity of Syrtis 
Major and Lunae Palus and are approximately 0.5° with respect to 
a smooth sphere. It may be significant that these slopes occur in 
conjunction with relatively strong echoes, and also in approximate 
coincidence with dark areas . If one draws a parallel with the 
distribution of snow in wintry regions of the earth (ignoring the 
temperature variation with altitude), it is neither the flat up- 
lands nor the flat lowlands that are barren, it is the steep sides 
of mountains. Thus, G, H. Pettengill (private communication) has 
suggested that it may be the relative absence of dust that causes a 
higher radar reflection efficiency as well as the darker visual 
appearance . 

It should be remembered in this discussion that both the cross 
section and the height represent an average over a region of about 
100 km in radius (1.7° are on the surface). This estimate of the 
size of the effective scattering region has been established from 
spectral observations by earlier workers and is confirmed by spectral 
data obtained in the measurements reported above as well. Thus the 
mean heights and slopes determined by Pettengill are large-scale 
averages and may hide substantially larger local variations. 

In comparison with the earth the sample of Mars at 21°N that 
has been studied during the opposition of 1967, displays far greater 
height dispersion. If the terrestrial ocean floors are included, 
this is no longer true. However, since Mars has only half the radius 
of the earth, the fractional large-scale height variation remains 
considerably larger. With only slightly more than one-third the 
surface gravity of earth and with considerably less erosion, Mars' 
greater relative surface relief may be understood. 
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RADAR CONTRIBUTIONS TO PLANETARY ATMOSPHERES 


Venus 

Cross section vs wavelength .- The suggestion has been made 
(see previous section) that the variation of the radar cross 
section of Venus with wavelength (Fig, 1) can be accounted for 
in terms of atmospheric absorption. Since the largest part of 
the echo is returned from a region around the subradar point 
where the rays penetrate the atmosphere nearly vertically we may 
neglect the curvature of the surface and to a first approximation 
write 

a = a Q exp(-2t) (1) 

where a Q is the intrinsic cross section (observable in the 
absence of an atmosphere) and T the "optical depth" defined in 

t = Jk dh (2) 

where K is the absorption coefficient per vertical height element 
dh. Evans et al. (1966b) suggested that the wavelength 
dependence of the absorption coefficient might be of the form 
K “ 1/X leading to a law of the form 

a = a Q exp (-a/ A) (3) 

and B. Lax in a companion paper (see Appendix Hi) has examined 
a law of the type K <= 1/A 2 leading to 

a = o Q exp(-b/A 2 ) . (4) 

These two laws are compared in Fig, 9, This figure differs 
from Fig. 1 in that in order to test these two laws (Equation 3 
and Equation 4) to the fullest we have reduced the error bars 
associated with the A = 12.5 cm measurement to the +1% vaiue 
given by Carpenter (1966) in his most recent paper. The error 
associated with the 23 cm point has been reduced from +2 db to 
to +1 db on the basis of a subsequent radar calibration using 
the Lincoln Calibration Sphere (a perfectly round machined metal 
sphere placed in earth orbit as a standard radar target) # This 
later calibration confirmed the earlier cross section measurement 
but reduced its associated uncertainty. The 70 cm point has been 
replaced by a revised value reported by Pettengill et al. (1967), 
viz. 14 + 7%. 


D-2 3 



CROSS SECTION (%) 


24 


20 

18 



10 100 1000 
X(cm) 


F igur e 9 The radar cross section of Venus and the curves pre- 
dicted by two simple models described in the text. The error 
bars shown in this plot have been refined over those given in 
Fig. 1 by the use of more recent values. 
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The curves drawn in Fig. 9 have been forced to pass through 
the value a = 1.2% reported by Evans et al . (1966) at 
X = 3.8 cm. By allowing a change in the value of a at this 
wavelength, yet remaining within the error bars shown, either 
curve could be made to pass through the error bars of the 
X = 12.5 cm point. 

In sum the present cross section measurements are incapable 
of distinguishing between the two most elementary wavelength 
dependencies for the absorption coefficient. It is possible 
that the wavelength dependence of the absorption is compounded 
with variations in the intrinsic cross section a Q with wave- 
length due to different penetration depths of the signals. Little 
can be said about this possibility other than to point out that 
in the case of the moon (where it is known to occur) the in- 
trinsic cross section is lower at all wavelengths. This suggests, 
as noted before, that Venus has a more compact surface and there- 
fore that this is not an important effect. Such a conclusion is, 
however, little more than surmise at present. 

Atmospheric scattering .- It is possible to argue that the 
low value for the cross section observed at 3.8 cm arises because 
the signals are no longer being scattered from the surface but 
from particles (e.g. rain or ice) in the atmosphere. Four 
arguments can be advanced against this hypothesis: 

1. The absolute doppler shift agrees with that predicted for 
the planetary surface to the limit of the experimental accuracy. 
Thus the particles would have to have equal upward and downward 
velocities with respect to the surface (or none at all) to be 

seen at the same velocity as the planet. 

2. The overall doppler width of the signals (Fig. 10) is 
precisely that computed from the known size and rotation rate 
of the planet. No signals are seen at frequencies outside this 
band that could be attributed to particles in turbulent motion. 
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Figure 10. A spectrum of echo power for Venus at 3.8 cm wavelength. Note that the 
peak echo occurs at the predicted doppler shift and all the power is contained between 
the two arrows which indicate the computed positions of the limbs (Evans, et al. 1966b) 



3. The features are seen in the spectra (e.g. Fig. 10) and 
these move in position within the spectrum in a systematic manner 
consistent with their being attached to the planetary surface. These 
features appear to be the same as those seen at longer wavelengths 
( see Fig . 5 ) . 

4-. Short pulse range measurements have been made essentially 
simultaneously (Evans et al. 1966b) at wavelengths of 3.8 and 
23 cm (Fig. 11). The results show that the nearest reflection 
point for the two frequencies is at the same range to within 
the experimental accuracy (+1 km) . 

In conclusion, if any scattering does occur in the atmosphere 
of Venus, existing experiments have failed to detect it. A 
radar operating at a frequency of 1 cm or less would be better 
suited to such a task because a) the atmospheric absorption to 
the surface would presumably be higher thus reducing the background 
"glare" of the planet, and b) the particle cross section 
would be increased by 200 times assuming that the particle radius 
is small compared with the wavelength. 

Limb darkening . * If the reduced cross section at 3.8 cm is 
to be accounted for in terms of atmospheric absorption, we would 
expect the absorption to be most severe near the limbs of the planet 
where the rays must penetrate the atmosphere obliquely. This 
would serve to steepen the angular scattering law P (4>) 3 g over 
that which would be seen in the absence of an atmosphere. If 
the one-way attenuation were A db , this modification can be ex- 
pressed in: 

10 log. _ P ( 4> ) - 10 log n F (<J>) = 2 A sec <f> db' (5) 

Thus, given some independent means of measuring F (<f>), it would 
be possible to deduce A directly, it should be noted that Eq. (5) 

is based on the premise of rectilinear propagation, thin atmosphere 
(i.e. scale height << planetary radius) and uniform absorption 
over the disc. A number of these assumptions can be questioned 
and in particular the last. If the poles are colder than the 
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Figure 11 . The echo power vs delay observed for the initial part of 
the return from Venus at two wavelengths. The delay to the reflecting 
surface is the same at both wavelengths to the limits of the experi- 
mental accuracy (+_ 1 km) . Further the shapes of the echoes are 
nearly identical. These results show beyond doubt that it is the 
surface and not the atmosphere that gives rise to the 3.8 cm echo 
(Evans, et al. . 1966b). 
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equator, then the total absorption there might be different. For 
example, for a dense CC^ - Ng atmosphere (Barrett 1961) in which 
the absorption is produced by collision-induced dipole moments 
(Barrett and Staelin 1964): 

A a (1/T 4 * 6 ) db (6) 

Accepting the limitations of this model, Evans et al. (1966b) 
searched for a limb darkening effect. In the absence of any way of 
measuring P (<|>)(i.e., the scattering behavior with the atmosphere 
removed) they were obliged to adopt as a model for P (<p) the curve 
observed at some longer wavelength, and then assume that this did 
not differ from the one that would be obtained at 3.8 cm. In the 
case of the moon such an assumption would be invalid. However, no 
systematic wavelength dependence in F (<j)) for Venus has yet been 
found (see preceding section). Two experiments were tried. In one, 
an echo power vs delay curve was constructed using pulses of various 
lengths, and is compared with the 23 cm results (Fig. 2) in Fig. 12. 
In the second, a mean echo power spectrum was constructed by 
normalizing the four best frequency spectra to a common relative 
scale (Fig. 13) and this was then compared with a similar mean 
spectrum obtained at 12.5 cm at J.P.L. (Muhleman 1965) in Fig. 14. 

In both cases the effects of limb darkening are evident. 

From the pulse measurement (Fig. 12) a value of 
2A = 5.5 + 2 db was obtained, and from the spectrum (Fig. 14) 

2A = 3.5 +_ 2 db. Thus in neither case was the limb darkening 
observed consistent with the amount required (2A 'v 10 db) to 
reduce the cross section to the low value observed. This discrep- 
ancy can result from a number of causes such as: a) the model 

may be oversimplified (as noted above), b) the intrinsic cross 
section a may be lower at 3.8 cm, and c) the angular scattering 
law P (tf)) may change with wavelength for Venus as it does for 
the Moon. This last eventuality is capable of introducing the 
entire discrepancy. In all probability, however, the true 
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Figure 12 . Comparison of echo power vs delay at 3.8 and 23 cm wave- 
length. The divergence at delays >1 m sec is thought to indicate 
limb darkening caused by the atmosphere of Venus (Evans, et al. 1966b) 
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Figure 13. The mean echo power spectrum derived by Evans, et al 
TltFFET'Tor 3,8 cm. 
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Figure 14 . Comparison of the curve of Fig. 13 (labelled X-band) 
with that obtained by Muhleman (1965) at 12.5 cm (S-band). The 
divergence is again thought to indicate limb darkening (Evans, 
et al. 1966b) . 
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explanation probably lies in a combination of all three effects. 

In conclusion this experiment would indicate that the one- 
way zenithal absorption of 3.8 wavelength signals is probably not 
less than 2 db while the cross section measurements (see above) 
imply that it cannot be more than 6 db. Refinements in these 
conclusions can be expected from further observations at 3.8 cm 
wavelength using the improved Haystack radar. Very desirable also 
would be new measurements at say X - 1,2, and 6 cm. 

Mars 

Radar experiments on Mars have, as yet, yielded no information 
concerning its atmosphere. 


DISCUSSION 

Surface Material of Venus 

The long-wave radar measurements suggest a dielectric con- 
stant of the order of 5 - 6 (see above). The value of e = 3.5 
reported by Carpenter (1966) should be regarded as suspect in view 
of the possible atmospheric attenuation at this wavelength (Fig. 9). 
The value e - 2.2 reported by Clark and Kuzmin (1965) from radio 
interferometer polarization measurements at 10.6 cm is likewise 
suspect. The curves of Fig. 9 imply the existence of between 
1.25 and 2.25 db zenithal absorption (one-way) at this wavelength. 
Thus at a position on the disc where the emission would appear 
polarized (<j> 60°) the one-way absorption will be between 2.5 and 

4.5 db. That is at least half the emission will be from the 
atmosphere and will be unpolarized. 

If correct, this argument serves to resolve the discrepancy 
between the radar and radiometric values for the dielectric constant, 
and incidentally argues against the cloud-greenhouse model of Sagan 
and Pollack (1965, 1967) which predicts essentially no absorption 
of 10.6 cm wavelength radio signals. This conclusion was reached 
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earlier by Hansen and Matsuchima (1967) who show rigorously that a 
surface dielectric constant of e = 6 is compatible with the fringe 
visibility observed by Clark and Kuzmin (1965) given an atmospheric 
attenuation as large as the values suggested above. 

Surface Material of Mars 

Strongest reflections from Mars are apparently associated 
with sloping terrain, and such terrain tends to be dark, though 
here the correspondence is less than perfect (Fig, 8). 

G. H. Pettengill (private communication) has suggested a model to 
account for this behavior in which it is supposed that the dark 
regions having high radar reflectivity are the sloping surfaces 
which are less densely covered in eroded material or debris. On 
the other hand, B. Murray (private communication) has suggested 
that the strong radar returns are from subsurface water or ice 
tables. This might then account for the extremely smooth appearance 
of these regions and at the same time account for an albedo dif- 
ference — the soil above being damp. Apart from the obvious com- 
ment that such regions might be expected to occupy the lowland 
ground rather than the slopes little more can be said, and it 
seems profitless to speculate further until more data is available. 

Atmosphere of Venus 

We believe that the atmosphere of Venus gives rise to between 

2 and 6 db atmospheric attenuation one-way at the zenith. Of the 
eight models for the atmosphere of Venus proposed by Barrett and 
Staelin (1964) to account for the microwave emission spectrum only 

3 can give rise to this much absorption. These are: a) the 

C0 o - N 0 model with a pressure of several hundred atmospheres, 

*■ 3 

b) the water vapor cloud model — 6 km thick with 1 gm/meter 

water vapor in equilibrium with liquid water at the cloud bottom, 

and c) the lossy (and therefore absorbing) dust model. 
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Sagan and Pollack (1965, 1967) favor a cloud-greenhouse model 
for the atmosphere in which CO^ largely provides the infrared 
opacity and is enhanced by cloud particles , and the latter provide 
the microwave opacity. This may be regarded as a combination of 
models a) and c) and has been criticized by Hansen and Matsushima 
(1967) on the grounds that the radar reflectivity and differential 
polarization of 10.6 cm emission (see above) indicate that the 
microwave opacity is concentrated most heavily near the planetary 
surface; and not in a low temperature cloud region. As an 
alternative they propose a variant of model c) in which the 

heat input is solely that of the internal heat of the planet. 

2 

The amount of dust required in the atmosphere above each cm 
of surface is then of the order of 10 gm. Despite these some- 
what bizarre propositions, the model is capable of matching 
both the microwave emission and radar reflectivity data. 

A major obstacle to this last model is the variability in 
the 3.8 cm cross section of Venus observed by Evans et al. (1966b) 
who suggested that there might be "clearings" in the atmosphere, 
i.e. regions in which the abundance of the microwave absorbing 
agent is substantially lower than the average. Such clearings 
would destroy the effectiveness of the dust blanket and would 
in any event be difficult to generate in an atmosphere where 
surface temperature differences (and therefore winds) mus't be 
largly absent. 

In the companion paper we examine the evidence for variations 
in the radar cross section of Venus and conclude that, in large 
measure, these can be accounted for by variations in the height and 
roughness of the terrain at the subradar point. 

FUTURE WORK 
Venus 

Important new results can be expected from a continuation 
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of existing programs — in particular the 12.6 cm J.P.L. observa- 
tions (hopefully with the 210 ft, antenna) and the 3.8 cm 
measurements at Haystack. Yet measurements at new wavelengths 
almost anywhere in the range £ 1 - 10 cm are urgently required. 
These offer the promise of determining the absorption coefficient 
vs radio wavelength for the Venus atmosphere and consequently 
will constrain the number of models that can be invoked to 
account for the microwave emission spectrum. In addition, ob- 
servations at the shortest of these wavelengths could yield 
direct evidence of particulate matter in the atmosphere. Accurate 
observations at long wavelengths (A ■> 2 m) are required in order 
to provide a reliable value for the intrinsic reflection 
coefficient of the surface. 

The new measurements at 70 cm, 12.6 cm and 3.8 cm promise 
to provide additional information on the location and nature of 
the anamolously bright scattering regions. Plans to employ a 
radar interferometer at 3.8 cm hold the possibility of providing 
a crude though useful radar map. 

Mars 

The extension of the height contour measurements (Fig. 8) to 
other regions promises to yield the first contour map for Mars. 

The variation of cross section with surface height and the 
associations of these two with surface markings are exciting 
results. With added data the meaning of these associations may 
be made clear. 

Long wave measurements (A ^ 2 m) would be extremely useful 
in conjunction with the 3.8 cm measurements. They might show 
for instance that the longer wavelength penetrates to a greater 
depth and encounters ice or water (as proposed by Murray above 
and Smoluchowski in Appendix J ) . 
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REFLECTION AND DIFFRACTION IN A BISTATIC 
RADAR OCCULTATION EXPERIMENT 

ABSTRACT 

The strength of the field reflected obliquely from a planet 
in an occultation experiment is computed for a number of para- 
meters, It is shown that a study of the reflected signal • component 
can provide information about the electrical properties of the 
planetary surface. Possible distortions of the diffraction 
pattern caused by surface roughness are also investigated 
and it is found that the likelihood of occurrence of surface 
undulations large enough to be observable is very small. 


-ii- 



TABLE OF CONTENTS 


Title Page 

Introduction E-l 

Evaluation of Reflection from the Surface. ........ E- 2 

The Diffracted Field .......... . ........ E - 7 

Conclusions E-10 

Reference ...... E-ll 

Appendix , 12 

LIST OF ILLUSTRATIONS 

Figure Title Page 

1 Spacecraft in Relation to the Planet ,. . . . . E-3 

2a Relative Reflected Field due to 

Reflection from Dielectric Planet as a 
Function of Position when Spacecraft is 

at Height 1.0 *a Above Surface E-5 

2b Relative Reflected Field due to 

Reflection from Dielectric Planet as a 
Function of Position when Spacecraft is at a 
Height 0.3 ‘a Above Surface ........ . . E-6 

3 Definition of Quantities Required to 

Compute Diffracted Fields E-7 

4 Effect of Uneven Edge on Diffraction 

Pattern , . . E-9 


-iii- 



REFLECTION AND DIFFRACTION IN A BISTATIC 
RADAR OCCULTATION EXPERIMENT 
T. Hag for s 

Lincoln Laboratory* , Massachusetts Institute of Technology 

July, 1967 

INTRODUCTION 

G. Fjeldbo and V. R. Eshleman (1965) have described methods 
whereby atmospheric properties can be derived by careful studies 
of phase and frequency of a signal received by a spacecraft being 
occulted by a planet. The method has been applied to studies of 
the atmosphere of Mars with considerable success. 

In that work little consideration was given to the possi- 
bility of deriving information about surface properties . The 
surface properties are indeed not of great importance in 
determining the signal properties as long as the spacecraft is at 
a large distance from the planet at the time of occultation. 
However, future occultation experiments both on Mars and Venus 
will employ spacecrafts which approach the planetary surface 
quite closely, and it is felt that an analysis of the amplitude 
pattern of the signal near occultation could well supply valuable 
information concerning surface material and possibly structure. 

To evaluate these possibilities some numerical calculations 
are carried out in this report in order to evaluate the relative 

* Operated with support from the National Aeronautics and Space 
Administration under Contract No. NSR-22-009-106 . 



strength of a direct and a reflected signal when the planet is 
smooth. Little progress has so far been made in evaluating the 
effect of surface roughness on the reflected signal but it is 
believed that the total strength is not strongly dependent on 
the roughness. With substantial surface undulations the 
occultation diffraction pattern will be affected. The conditions 
for this to occur are also worked out in some detail. 

In the process of this work, inspired by several good sug- 
gestions from D. L. Anderson, California Institute of Technology, 
an exact method was developed to determine the refractivity 
profile of an atmosphere from the Doppler-offset versus time 
behavior of the signal. The method apparently is well known in 
seismic sounding techniques . 

EVALUATION OF REFLECTION FROM THE SURFACE 


The geometry of the spacecraft in relation to the occulting 
planet is shown in Fig. 1. The spacecraft is assumed to move in 
the xy plane and to have coordinates R, as shown. The direction 
to the observer (at infinity) is taken to be along the negative 
y-axis. We shall assume that the spacecraft is transmitting. 

The signal at the observer can be considered to consist of two 
parts, one derivable from the free space field on the plane 
y = 0 exterior to the planet and another derived from the field 
distribution set up on the planetary surface. When the distance 
R is much larger than the planetary radius a the former is 
dominant , but for smaller R the reflected component may become 
appreciable. In order to assess the importance of this component 
we compute the reflected signal strength by Huygens' principle: 
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where 

E^R-j^ = free space field strength at distance R^ 

d£ = surface element = a 2 d^ sin0 d0 
n = surface normal 

R(cosa) = reflection coefficient at an angle of incidence a 
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R = distance from center of planet to observer 


The stationary phase point ip_ 0 is found by equating the partial 
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derivatives of (R^ + R 2 ) with respect to ip and 0 to zero. It is 
found that: 

0 = 90° and §• cosip = cos (0 - 2ip ) (2) 
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Expression (3) was evaluated on a digital computer for <j> ranging 
from -90° to <j> s Arc cos (a/R) at which point geometric 
occultation occurs. The calculations were carried out for a 
dielectric planet with dielectric constants of 4.0 and 6.0, and 
for values (R/a) of 1.1, 1.3, 1.5, 1.8 and 2.0, Some of the 
results of the calculations are plotted in Figs. 2a and 2b. As 
can be seen the reflected signal may well be appreciable. The 
very different behavior of the two principal linear polarizations 
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Fig. 2a Relative Reflected Field due to Reflection from 
Dielectric Planet as a Function of Position when Spacecraft 
is at Height 1.0 'a Above Surface. 
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Fig. 2b Relative Reflected Fie 
Dielectric Planet as a Function 
is at a Height 0.3*a Above Surf 
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should also be observed. Given enough signal strength a close 
study of the modulation as well as the polarization of the 
received signal hence can provide valuable information on the di- 
electric properties of the planetary surface. 

THE DIFFRACTED FIELD 

In addition to the reflected field evaluated in the preceeding 
section there is a diffracted field which can be determined 
approximately from Huygens ' principle by integrating the free- 
space fields set up by the transmitter over the plane y = 0 
exterior to the planetary body, see Fig. 3. 



Fig. 3 Definition of Quantities Required to Compute Diffracted 
Fields . 


E-7 


The free-space field over the plane y = 0 is: 


0 O <X, z) = 




* — Y 

cos<j>) +R 


sin <|> + z 
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(4) 


As long as a » /X R sin<{> we may regard the limb of the planet 
as a straight edge normal to the xy plane at y = 0 , x = R cos<j>-a. 
In terms of the Fresnel integrals one finds for the field at the 
observer: 

E -ikCR +R sin<$>) C, 

^o = 2l~ |f(l-i)+C(A) - iS (A )J (5) 

where : 

A = /2~ (R cos<t>-a)/v'RX sin<(> 
and where the Fresnel integrals are defined by: 

C(x) = J* cos ^ S(x) = J sin ^ <p 2 Jd<t> 

Expression (5) which is well known from previous work is to be 
combined with (3) to give the total received field at the 
observer. 

We next investigate the possibility that surface roughness 
may play a significant part in determining the properties of the 
diffracted field. To simplify matters we consider the edge to 
be straight in the mean and to have random deviations super- 
imposed according to: 

x = a + f(z) (6) 


where f(z) is a zero-mean random function of z, see Fig. 4. 
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Fig . 4- Effect of Uneven Edge on Diffraction Pattern 


Using the same approximations and notations as before we obtain: 
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where now; 


ACz) = /2 |R cos<J> - a - f(z)j =A-/Tf(z) 

Expanding the Fresnel integrals about A(z) = A for reasonably 
small values of f(z) (see Appendix) one obtains from (7): 
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Little further progress can be made here unless a specific 
functional dependence f(z) is assumed. However, a crude general 
conclusion can be drawn. Appreciable corrections to the straight 
edge result can only appear when f (z) is of the same order of 
magnitude as /X R sin<Js, i.e. the extent of the central Fresnel 
zone, and this variation of f (z) must occur in a distance of 
order ^ /A R sin<J> . Numerical estimates indicate that even at 
the shortest wavelengths such variations are quite unlikely to 
occur. Therefore, information about surface roughness does 
not seem to be readily available from studies of distortions 
of the diffraction pattern. 

CONCLUSIONS 

The present note has shown that the signal component which 
is reflected from the planetary surface may be appreciable and 
that it can be used to derive information about planetary 
surface properties of considerable interest. 

The note has also shown that distortions in the diffracted 
field caused by surface roughness are very unlikely to be 
observable under conditions expected to prevail on planetary 
surfaces . 
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APPENDIX 


NOTE ON AN APPROXIMATION RELATING TO FRESNEL INTEGRALS 


Let us define the function F(x) by: 

x 


F(x) = C(x) - i S(x) = \ exp 

o 


= J (-i f * 2 j 


d4> 


( 9 ) 


We would like to derive a simple approximation to F(x+Ax) in 
terms of F(x) and elementary functions of Ax. To do this we 
first write FCx) as: 
x 2 


(x) = J § exp {- i j tj dt/Zt = 


x 


= | | dt exp i | t - | log tj 


( 10 ) 


Hence 


(x+Ax) 


= | e xp ^-i | 
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Ax 2 

I Jk 

I , \ 

, 2 
log X 

J \ dt ' exp 

-t* (i \ + 1/2 X 2 ] 

J 

o 

l 1 7 


^ f + 1/2 x 2 )^ 


(1/X) + i TT X 


( 11 ) 


This expression is valid over a region considerably greater than 
that of a Taylor expansion of F(x) about the same point. 
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WATER VAPOR AND ICE IN THE MARTIAN ATMOSPHERE 


ABSTRACT 

The infrared measurements of water vapor content in the 
Martian atmosphere together with the atmospheric temperatures 
inferred from the Mariner IV occultation experiment seem to indi- 
cate that the atmosphere is saturated with water vapor. Since 
winds and clouds have been observed on Mars , it is reasonable to 
assume that there is considerable vertical mixing in the lower 
atmosphere. Moreover, the warmer Martian surface will act as a 
reservoir to maintain a saturated or slightly supersaturated 
water vapor content in the atmosphere above it. Thus the model 
for water vapor content is one in which the lower atmosphere is 
saturated with water up to a height of some tens of kilometers. 

In this region water ice crystals will form whose size is de- 
pendent on temperature and height. The existence of these ice 
crystals can perhaps explain the blue haze, diurnal clearings of 
haze, the white clouds, blue clearings, and the clouds observed 
on the limb in the Mariner photographs and in other earth-based 
photographs. Water ice crystals cannot be invoked to explain the 
yellow clouds or the polar caps. 
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WATER VAPOR AND ICE IN THE MARTIAN ATMOSPHERE 

H. Heffner 
July, 1967 


INTRODUCTION 

The existence of particulate material in the Martian atmos- 
phere is attested to by a wide range of visible phenomena. They 
include the almost constant presence of "blue haze" and the 
formation, movement, and dissipation of several types of clouds. 
These clouds fall into classes termed: yellow clouds, white 

clouds, and blue clouds depending upon the optical wavelength 
where greatest contrast is observed. The blue haze refers to 
the obscuration of surface detail when the planet is photographed 
over a wavelength range in the blue and ultraviolet. 

Many explanations have been tendered to account for the 
various types of clouds and for the blue haze. They have 
included hypotheses of great dust storms , of particles of water 
ice, and of particles of solid CC^- Each of these hypotheses was 
put forth before the 1965 Mariner IV occultation experiment 
revealed the pressure and temperature profile of the Martian 
atmosphere, and little re-examination of the models in the light 
of this new information has been carried out. 

As a result of the occultation experiment''' , we now believe 
the Martian atmosphere is almost wholly composed of CO2 with a 
surface pressure 4.1 to 5.7 mb. at the immersion point and a 
scale height of 8.5 to 9 km. The experiment seems to indicate 
an isothermal atmosphere to some 20 km or so above the 
immersion point having a temperature of 160°K (Fjeldbo , private 
communication). Since the surface temperature is some 60° to 90°K 
higher than this , we expect a thin atmospheric boundary layer 
where the temperature varies rapidly from perhaps 250° K at the 
surface to 160°K and then remains essentially constant with 



height for perhaps the next 20 km. At the point of immersion 
(local summer night time) the surface pressure was approximately 
9 mb with a scale height of 10.5 to 14.2 km. The atmospheric 
temperature was estimated to be 250°K to perhaps as low as 220°K. 

Observations of water vapor absorption lines in the Martian 
atmosphere are rendered difficult by the Earth's own atmosphere, 

ii 2 

however, the experiment of Kaplan, Munch, and Spinrad was 

sensitive enough to yield an estimate of the precipitable water 

above the Martian surface of (14 ± 7) x 10 gm cm . If one 

assumed a constant mixing ratio in the Martian atmosphere, then 

the scale height of 9 km would imply a surface density of H„0 

13 - 3 ■ 

molecules of (2.75 +.1.37) x 10 cm and for a temperature of 
175°K, a surface pressure of 0.663 +_ 0.332 microbars. 

At this temperature, the saturation vapor pressure for ice 
is about .08 microbars, or almost an order of magnitude less 
than the pressure derived. Johnson has cited this discrepancy 
to question the accuracy of the measurement. It is reasonable 
to believe, however, that the summer hemisphere exhibits higher 
temperatures so that an integrated disc average would yield the 
amount of water vapor measured. 

The value is sufficiently close to the saturated value that 
it suggests further investigation of the implications of a model 
Martian atmosphere which is saturated with water vapor and which 
contains ice crystals in a variety of sizes and numbers to account 
for most of the clouds and haze observed. 

THE MODEL FOR ATMOSPHERIC WATER VAPOR 

If our experience of Earth conditions is any judge, the 
outgassing of the Martian surface must have released immense 
amounts of water, yet very small amounts are observed in the 
atmosphere. The explanation probably lies not only in the fact 
that great amounts of water vapor probably escaped, but because 
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of the low temperature of the planet a considerable amount of 
ice exists on and below the surface. Even though this ice is 
likely to be in capillary fissures or covered with sandy material, 
it can act as a reservoir for water vapor to be given off to the 
even colder Martian atmosphere. Because the temperature declines 
rapidly with height in the narrow boundary layer just above the 
surface, the vertical mixing action of wind will attempt to keep 
the water vapor pressure above the saturation value for some 
distance up. The atmosphere of course cannot remain appreciably 
supersaturated for long until ice crystals begin to form and to 
grow in sufficient size and number that the partial pressure of 
H2O is reduced to nearly the saturated value. As a result, we 
expect the partial pressure of water vapor to stay essentially 
constant at approximately the saturation value for a region 
extending perhaps to 25 km above the boundary layer. Mixed in 
will be small ice crystals of sufficient size and number to keep 
the partial pressure of HgO near the saturation pressure in the 
face of vertical mixing from the boundary layer with its 
relatively high water vapor content. The value of 25 km is taken 
for the extent of this region principally because the Mariner IV 
data show this to be the region of constant temperature , More- 
over, Leighton has interpreted a light area on the limb in a 
Mariner IV photograph as being a cloud with its top at 25 km. 

It should be noted that even at saturation, water vapor 

would make up an almost negligibly small fraction of the 

total atmosphere. With the above model, even at 25 km the 

density ratio of water vapor to C0 0 is less than 10” ^ at 160°K 

-7 

and at the surface is less than 10 

PARTICLE SIZES IN THE MARTIAN ATMOSPHERE 

From extensive Russian observations of limb darkening, 

”5 

Opik quotes the following table of optical characteristics of 
the Martian atmosphere and surface. 
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TABLE 1 

TRANSMISSION COEFFICIENT (p) AND SPECIFIC SCATTERING OF THE 
MARTIAN ATMOSPHERE (a) AND REFLECTIVITY OF THE SURFACE (s). 


Wavelength (A°) g a s 


4600 (Blue) 0 

. 33 

0.20 

0.25 

5200 (Green) 

.54 

.22 

.29 

5430 (Green-Yellow) 

. 60 

.23 

. 34 

5800 (Yellow) 

.69 

.24 

.40 

6400 (Red) 

.74 

. 20 

.53 

6 

Brooks plots the curve 

shown 

in Fig. 1 for the 

total albedo 


geometric reflectivity of Mars showing a minimum in the vicinity 
of 3500 A° . 

From these data we can gain some rough ideas concerning the 
size of the particles and their density. The intensity in a 
scattering atmosphere goes as 


I = 


-tta Q ext * 


( 1 ) 


where N is the volume density of the particles, a their radius 
and Q ext is the coefficient of extinction. Figure 2, taken from 
van de Hulst^ shows computed values of Q ext for spheres with 
index of refraction of 1.29. The curve peaks for k near zero at 
a value of 2ira/A - 6.5. Taking this point to be the minimum 
reflectance point of Fig. 1 at 3500 A°, we find a particle 
diameter of 2a = 0.7 microns. By matching the attenuation to 

tt n 7 

the value given by Opik for 4600 A we obtain Ni = 2.2 x 10 

particles per square centimeter. This number is lower by an 

0 

order of magnitude than that of Greenspan and Owen who derive a 
particle density required to bring into agreement pressures 
derived from both visible and ultraviolet measurements. 

Assuming the particles are ice with diameter of .7 microns, we 
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GEOMETRIC REFLECTIVITY (p) 



2000 3000 4000 5000 6000 7000 

WAVELENGTH (X, ANGSTROMS) 


Figure 1 The geometric reflectivity of Mars . Solid line, 
relative data adjusted to 0.04 at 3400 A; +, absolute 
reflectivity determined by comparison with $ Leo, plotted 
independently. The dashed lines below 3400 ^ represent the 
error range applied to the relative data. 
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find their total mass to be 4 x 10 - gm/cm . This value is 
less than one hundredth of the observed amount of precipitable 
water vapor. 

Ice crystals of this small size will stay suspended 
indefinitely. Being much smaller than a mean free path, their 
fall rate is governed not by Stoke’ s Law but rather by non- 
viscous flow considerations. The terminal velocity for these 
particles of radius a smaller than a mean free path is given by 


v 


= lR.S a — 

3m 'nv 

Cl 


( 2 ) 


where p is the density of the particle, g is the acceleration 
due to gravity and m* , n and v are respectively the mass, 
number density, and average thermal velocity of the surrounding 
gas molecules. For particles of 0.7 microns diameter in the 
Martian atmosphere, the fall rate is approximately 0.04 cm/sec 
and is thus negligible in comparison to effects of winds. 

The growth of the ice crystals is governed by diffusion. 
Thus, the rate of growth is given by 

da 1 dm _ ^Pw - ^i^ 

dt “ dt ’ PK„T a 

HTTa 0 W cl 


where D is the diffusion constant for water vapor, p r and p. 

are the ambient water vapor pressure and the vapor pressure of 

the ice crystal respectively, p the density of ice, T the 

ambient temperature, and R the specific gas constant for water 

w 

vapor. 

The vapor pressure of the ice crystal will be determined by 
its temperature according to: 

p. s Ke" L/R w T i (4) 

x 

7 

where L is the heat of sublimation of water (2834 x 10 erg/gm) , 
T^ is the temperature of the ice crystal, and K is a constant of 
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proportionality. Numerically , this equation becomes 

P i = (3.63 x 10 7 ) 10“ 2670/T i (5) 

The temperature of the ice crystal is determined by the 
rate at which heat is transferred to or away from it. There are 
four main effects. First, heat is absorbed from the Sun at a 
rate which is determined by the solar constant on Mars, S the 
size of the crystal, and its absorption coefficient. If we 
consider a spherical crystal of radius a with an average 
coefficient of extinction in the visible due to absorption of 
Q a b s » then the rate of solar energy absorbed is 


df^ 

dt 




( 6 ) 


Heat is also radiated, primarily in the infrared and so there is 
a heat loss term of 


dH 2 

“dt 




(7) 


In general, we expect Q ' a ^ g to be smaller than Q since the 
particle size is smaller compared to the wavelength in the 
infrared than in the visible region. However, the emissivity 
of ice is much greater in the infrared than in the visible 
region, so that the net effect is probably one of increasing the 
value of Q’ abs over Q abs - 

Heat is also gained or lost through sublimation or evaporation. 
This heat flow rate is 


dH 3 _ , dm ‘" ,a 

dt “ U dt " R T 

w 


( 8 ) 


Finally, heat is lost through convective cooling by the pre- 
dominantly CO 2 molecules striking the crystal. Assuming the 
accommodation coefficient to be unity, this rate is 
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(9) 


dH 4 

dt"~ 




(T) v 

Cl 


(T) 



where (T) is the pressure of C0 2 at the ambient temperature 

and height of the particle. 

At equilibrium, we expect both ^ and ^ to be zero. Since 
cooling by thermal radiation is negligible compared to convective 
cooling, these two conditions yield and 


V T . AT 2 W 

"T 3 P co (t)v 4 (I) 


Using the vapor pressure-temperature relation, we find 


Pw 

Psat 


2L AT 
R T T 


w 


4LQ abs S/3 R„T p c02 <T)v a <T) 


( 11 ) 


Here p gat is the saturation vapor pressure of water vapor at the 
ambient temperature and p w is its actual ambient partial pressure 

The extinction coefficient due to absorption, Q for 

g B.DS ? 

small particles is proportional to the radius . For ice crystals 
such that 27Ta/X is small 

W “ * Im (m - 1) (12) 


where m is the complex dielectric constant. Thus, if one knows 
the spectral variation of absorption of ice and the degree of 
supersaturation, (p /P ,), then he. can compute the equilibrium 
radius of the ice crystals which are formed. 


Several qualitative inferences can be drawn from the above 
relations. First, the crystal size should decrease with height. 
This effect results from the slower cooling due to reduced C0 2 
pressure and the correspondingly smaller radius which the 
particle reaches before its temperature rises sufficiently for 
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its vapor pressure to equal the ambient value. Further, we 
expect that as the atmospheric temperature rises, the degree 
of super saturation decreases and so does the particle size. 
Conversely, a decrease of atmospheric temperature causes crystal 
sizes to increase. 


ICE, HAZE AND CLOUDS 

This model of the Martian atmosphere is appealing since it 
can explain many of the visible features of Mars, As has been 
mentioned, the attenuation of the atmosphere can be attributed 
to the almost continual presence of minute ice crystal less 
than one micron in diameter. This blue haze is occasionally 
observed to clear over large areas in a matter of a few hours 
only to reform after a period of hours or days. This phenomenon 
seems impossible to explain if one assumes the scattering particles 
are dust and that clearings simply represent a settling of the 
windborne dust on the surface. Settling times for particles 
small enough to have the proper scattering characteristics are 
far too long to allow such rapid clearings. The model which 
accounts for the blue haze by ice crystals can allow for rapid 
clearings simply by atmospheric temperature increases which act 
to reduce the crystal size or conceivably, to cause them to 
disappear completely. 

In addition to the blue haze there is a generally whitish 
haze which varies diurnally, tending to clear off by Martian noon 
and to reform in late afternoon^ , Such a characteristic is pre- 
cisely what one would expect were somewhat larger ice crystals 
responsible. Relatively large crystals would form in the cold 
atmosphere toward late afternoon and dissipate as the sun heats 
the atmosphere each noon. 

Additional evidence tending to confirm the water vapor 
model arises from the polarization measurements of the white 
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clouds made by Dollfus'^' 1 '. From his studies, Dollfus identified 
these clouds as being ice crystal mists. 

Water ice cannot explain either the polar caps or the yellow 
clouds. Because of the very small amount of water vapor present 
even at saturation, unreasonably high wind velocities would have 
to be invoked in order to account for the transfer of water from 
pole to pole as the seasons change * It is much more reasonable 
to believe that the polar caps are composed of solid CC^ and 
that the clouds observed immediately over them are CC >2 crystals 
rather than 1^0 . 

The yellow clouds must be composed of considerably larger 

particles than those responsible for the blue haze and white 

clouds. This fact can be inferred in part from the characteristic 

yellow color which would be difficult or impossible to obtain 

from Rayleigh or Mie scattering, and in part from the settling 

time determined by the cloud durations. The particles must be 

12 

10 microns or more in diameter and the color must be charac- 
teristic of their composition. The standard explanation suggests 
they are dust particles transported upward by winds, by volcanic 
activity or by meteoric impact. It appears to be impossible to 
raise sufficient dust by meteor impact to account for the extent 
of the clouds observed if one assumes a meteor size-frequency 
distribution compatible with the crater counts made on the 
Martian surface . One is then left with the hypothesis of 
volcanic origin — for which there is no evidence — or of wind. 

The latter seems more probable. 
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THE VENUS ATMOSPHERE: 

FUTURE EXPERIMENTS AND PRESENT KNOWLEDGE 

ABSTRACT 

A review of present knowledge of the parameters most 
important to an understanding of Venus stimulates vital questions 
about the systematics of the lower atmosphere and its interaction 
with the surface. An atmospheric sounder (or entry probe) is 
proposed as a logical initial investigation to help satisfy most 
of these questions . General requirements for sensors and other 
devices to sample various characteristics of the atmosphere 
during the probe are included. The paper is concluded with a 
series of recommendations concerning the composition, design 
development, and operation of the proposed entry probe that are 
hoped will be helpful to planners considering further investigations 
of the planet Venus. 
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THE VENUS ATMOSPHERE: FUTURE EXPERIMENTS AND PRESENT KNOWLEDGE 

J. J. Hopfield and R. A. Phinney 
July, 1967 


INTRODUCTION 

This paper begins with a brief review of the present 
knowledge of the parameters most important to an understanding 
of Venus. The principle thesis of this paper is that in view of 
this knowledge the most crucial questions about Venus concern the 
systematics of the lower atmosphere and its interactions with the 
surface. The answering of the question "what is the physical and 
chemical nature of the atmosphere, clouds, and surface of Venus" 
should motivate the planning of missions to Venus in the immediate 
future. Recommendations arising out of this paper are given in 
the discussion section. 

In planning missions to the neighborhood of Venus it is 
essential to emphasize the great importance of an atmospheric 
sounder relative to a fly-by or an orbiter. Admittedly, any 
Venus probe which uses the "bus" method will have available a 
relatively large spacecraft which does not enter the atmosphere. 

We do not quarrel with using the bus as an experimental platform, 
but feel that first priority should be given to making it as 
effective as possible in support of the drop-sonde. We want 
to deliberately contrast this paper, with its strong emphasis on 
priorities, with other documents which regard all potential 
experiments as being equal in some sense, subject only to 
engineering constraints. (Reference JPL reports). In spirit, we 
follow the approach of the 1965 Space Science Board summer study, 
(Space Research: Directions for the Future; NAS-NRC publication 

1403; 1966), but hope in addition to add to the general priorities 
stated therein by specific consideration of the most crucial 
experiments . 



The experiments on a Venus mission can be divided into two 
broad classes. One is concerned with the physics of a tenuous 
atmosphere and ionosphere and its responses to the solar radiation 
and particle flux environments. The second is concerned with the 
lower atmosphere and its relation to the solid planet. There is 
extremely little which can be conclusively derived of fundamental 
interest to the second set of experiments from the first set. In 
addressing ourselves to the logical next step of this second 
problem, of interest to geology, astronomy, cosmology, and biology, 
we can virtually ignore the state of knowledge of the upper 
atmosphere . 

The only experiment likely to give further information on the 
basic atmospheric problem of Venus before the 1972 Mariner flight 
is the S-band occultation experiment on the current Venus probe. 
Even if successful, it cannot represent the major breakthrough 
for Venus that it did for Mars. While many detailed observations 
can be made on the ground, their nature has not and will not 
permit a decisive breakthrough on the Venus atmosphere problem. 

The first part of the following material discusses the state 
of knowledge of the physics and chemistry of the Venus atmosphere 
as an aid to understanding: 

a) why the atmospheric problem is not solvable from remote 
sensing ; 

b) why the atmospheric question is at present the important 
question on Venus; and 

c) what parameters are of relevance. 

The second part discusses the instrumentation of a sub-sonic 
atmospheric sounder designed to answer the most relevant 
questions . 


PRESENT KNOWLEDGE 

The gamut of conventional photographic , spectrometric , and 
radio and radar astronomy studies have been performed on Venus. 
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The most important conclusions of these studies are: 

1) Venus is covered by "clouds" at an altitude of about 

35 kilometers above the surface. This cloud cover has structure, 
but the surface of the planet has never been observed through 
rifts in the clouds. 

2) The cloud-top temperatures are of the order of 250°K. 

The cloud-top pressure is in the .05 - 1 atmosphere range. There 
are several independent methods of determining these "pressures" 
and "temperatures" from ground-based astronomy. (These methods 
do not directly measure temperature and pressure. Differences 
between the results obtained by different methods are expected 
and do occur.) There are observed (and probably real) temporal 
variations of the deduced temperatures and pressures. 

3) The combination of radar and radio astronomy show that 
the mean surface temperature of the planet is about 620°K, with 
a 150° day-night difference, and a relatively small north-south 
difference . 

4) CO2 is a sizable constituent of the atmosphere, although 
it is probably diluted by a major and yet unidentified component 
such as N 2 or A. 

These important questions arise from the observations: 

a) What is the physical, chemical, and geological nature 
of the atmosphere, clouds, and surface? 

b) Is this nature of the planetary surface a logical 
consequence of minor differences between the initial conditions 
of Earth and Venus, or were there important differences in the 
initial conditions at the time and place of the formation of 
Earth and of Venus? 

c) Has (or will) Venus go through a thermal and chemical 
configuration hospitable to the development of life? 

Much speculation has been written on these three questions. 

In spite of this, the questions are far from being satisfactorily 
answered. The logical dependencies of these questions necessitate 
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answering the first one first. The Venus space missions in the 
near future should be designed for the explicit purpose of 
answering this question. 

The Venus atmosphere keeps the planetary surface at a 
temperature almost twice that to be expected from an atmosphere- 
free planet the same distance from the sun. The origin of this 
"greenhouse" effect is a major puzzle. ^ 

Super-adiabatic lapse rates (vertical -temperature gradients) 
in an atmosphere are not stable. As a result, there is an 
upper limit on the mean temperature gradient in an atmosphere. 

For all likely gasses, this puts an upper limit of 7° - 10°/ 
kilometer on the Venus atmospheric temperature gradient. This 
lapse rate is consistent with estimates of the thickness of the 
atmosphere. It immediately leads to an estimate of the surface 
pressure of the order of 10-20 times the cloud-top pressure, 
i.e. ^ 1-20 atmospheres. The uncertainty in thickness and lapse 
rates lead to a large uncertainty in the surface pressure — a 
design range of 1 - 100 atmospheres should be allowed. A common 
feature of all models (in order to fit the observed data) is a 
thick, quasiadiabatic atmosphere. 

Cooling of the lower atmosphere (and surface) by radiation 
(and convection) tend to reduce the temperature gradient in the 
lower atmosphere. A high infrared opacity is essential to avoiding 
such radiative cooling. The three model mechanisms which have 


^"Throughout this paper we will presume the surface temperature 
is ^ 650°K. Although this is by far the most likely interpretation 
of the relevant data, there is in an indirect line of argument 
always some small possibility of gross misinterpretation. This 
possibility seems small enough that it need not be discussed here. 
On the other hand, it would be both scientific and engineering 
folly not to design experiments with sufficient dynamic range to 
cope with a possible lower surface temperature and much lower 
greenhouse effect. 
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been proposed for such opacity are: 

1) gaseous C0 2 and H 2 0 at reasonable (1-10 atm) pressures. 

2) Pure CO 2 with a high pressure atmosphere. 

3) Dense clouds of particulate matter having appropriate 
particle sizes and opacities. 

With a planetary surface temperature of 650° and a cold top 
of its upper atmosphere, general heat conduction processes will 
necessarily result in some upward heat flow by convection, 
radiative conduction, and/or precipitation and evaporation. For 
heat balance a source of heat at the surface is required. The 
models proposed for this heat source include 

1) penetration of solar radiation in the visible 

2) internal heat release 

3) frictional heating of the surface by high velocity winds 

4) a "thermal engine" driving adiabatic compression. 

The first of these is the mechanism by which the "greenhouse" is 
driven on the earth. The second is feasible if Venus is at 
present passing through a period of intense volcanism. The third 
and fourth note that the winds produced by insolation differences 
will drive circulation patterns which themselves can be used to 
heat the surface. The third mechanism uses frictional heating from 
these winds, while the fourth uses this circulation to adiabati- 
cally compress high altitude air (raising its temperature) while 
carrying it downward toward the surface. 

All models take into account some form of cloud cover. The 
three cloud materials most seriously suggested are C0 2 (solid), 

H 2 0 (solid or liquid), and non-volatile inorganic solid dust 
particles . The presumed role of the clouds ranges in models from 
1) transparent in IR; barely opaque in visible to 2) very dense 
and opaque at all wavelengths. The observed albedo of Venus is 
explicitly due to the cloud cover and determines the heat balance 
and "infrared temperatures" of the outer atmosphere. 
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The present atmosphere can be regarded to be in (perhaps 
evolving) equilibrium. On a short time scale, the atmosphere 
must by some mechanism maintain its heat balance, cloud cover, and 
a stable composition in spite of perturbations. (There is such a 
stabilization, for example, in the earth cloud cover.) The partial 
pressures of some of the constituents are also likely to be 
quantitatively understandable from the appropriate equilibrium. 

For example, the C^/CC^ ratio on the earth is maintained by plant 
life. The C0 2 pressure on Mars is probably regulated by the 
temperature of its north polar cap. On the high temperature 
Venusian surface, carbonate-silicate reactions are likely to 
buffer the C0 2 partial pressure. The relation between the 
surface, cloud, and gaseous atmospheric components is thus a 
problem any atmospheric model must try to understand. 

The variety and cleverness of ground-based observations of 
Venus have yielded an incredible amount of detailed information 
on Venus. The interpretations of this information are chiefly 
indirect. The great diversity of models now still "in agreement" 
with the observational data is a consequence of this indirect 
relation between theory and experiment. At this point, it is not 
likely that a few more years of conventional observations will 
produce any definitive progress in the understanding of the Venus 
atmosphere. When, on the other hand, a correct zero-order 
model for the atmosphere is given, presently available data will 
add much to the understanding of this atmosphere. 

In the period of time between now and the sending of the next 
interplanetary mission to Venus (Mariner 1972 or Voyager 1973) 
the following fundamentally new information is likely to be added 
to our knowledge of Venus. 

1) An S-band occultation experiment yielding the atmospheric 
number density and temperature profile in the upper atmosphere. 
Interpretational problems and atmospheric attenuation are unlikely 
to allow the extension of the altitude of this profile below an 
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altitude of 5 - 20 kilometers, or to pressures over 2-3 atmos- 
pheres . 

2) A map of the radar reflectivity with a resolution of the 
order l/40th of the planetary diameter will be made. 

3) The atmospheric attentuation of Venus at radar wavelengths 
will be known at several more wavelengths and with greater 
precision. 

In addition, much more information on the spectroscopy of the 
cloud top layer and its variability, the distribution of surface 
heights and surface slopes, the surface distribution of radio 
emission., and the physics of the Venus ionosphere will be 
available. Only the S-band occultation, however, makes any real 
progress toward a detailed understanding of the bottleneck in 
the scientific study of Venus, namely the understanding of the 
general atmospheric problem. 

ATMOSPHERIC SOUNDER 

For the sake of this discussion, we will posit an entry 
probe whose experiments weigh a few tens of pounds , and which is 
slowed to subsonic fall at an atmospheric pressure of about 5 
millibars. It is possible to design a group of experiments which 
will answer the present major questions about the atmosphere. 

The probe will require at least two minutes to pass to the 1 
atmosphere level, and several minutes more to reach the solid 
surface, depending on the drag coefficient of the probe and the 
atmospheric pressure profile. No special provision is made for 
these experiments to be usable at ionospheric or supersonic 
altitudes. It is assumed that the data will be telemetered 
essentially in real time, and that the probe will land "hard." 

(It would be desirable for the probe to survive landing as a 
working device, but this must be given low priority at this 
time.) The (probable) high surface temperatures present obstacles 
to the survival of conventional instrumentation; the short-term 
viability of the package at high temperatures is essential, 
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however, since the most important single objective is to determine 
the position, pressure, density, and temperature within a kilometer 
of an identifiable solid surface. 


Engineering Sensors 

We now discuss the foliowing engineering sensors which are 
required more or less independent of other experiments . 

1) Accelerometer . - To be integrated to give the position of 
the probe. If the dynamic response of the probe is well enough 
known, then an effective density can be inferred. With another 
instrument to measure density, the accelerometer data can be used 
to infer vertical winds. Only an x axial (z) component is needed 
for the important altitude determination*, the other two 'would 
give horizontal winds, and can probably be justified only for 
their diagnostic engineering value. 

2) Radar Altimeter . - This is needed to fix the relative 
position of the probe and the planetary surface. 

3) Communication Link. - It is possible to phase lock the 
carrier to an oscillator at the master station (bus or ground 
station) . Integration of the doppler shift gives an optical 
path distance to the probe. It is not clear that this data adds 
anything essential to the information from the accelerometer and 
altimeter, although redundance may be desired. There seems to be 
no real reason to want to measure refractivity of the atmosphere 
if the variables of real interest are to be determined directly. 

4) Thermal Sensors . - Thermal regulation of the probe is 

required both to cool the infrared sensors and to preserve the 
electronics at the high near-surface temperatures. The necessity 
of having engineering thermal data suggests that all thermal 
sensors be designed as a compatible system, including those 
sensors designated specifically to measure the local atmospheric 
temperature . 
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Thermodynamic Sensors 


We next describe the group of sensors for determining the 
local thermodynamic properties of the atmosphere. Little 
elaboration is needed; the instrumentation is well understood. 

1) Atmospheric Temperature. - One or more temperature 
sensors for determining the kinetic temperature of the atmosphere, 
appropriately shielded from radiation and appropriately placed so 
that the heating produced by the probe can be determined. 

2) Atmospheric Pressure . - One or more pressure sensors. By 
situating one at a point where the relative motion of probe and 
atmosphere are small, and one directly in the airstream for a 
dynamic measurement, both pressure and airspeed would be obtained. 

3) Atmospheric Density, - Direct density measurement. This 
can be done easily by measuring the beta absorption of a short 
column of "■air". In conjunction with P and T, this gives a mean 
molecular weight (M) . 

4) Velocity of Sound. - A measurement of the acoustic 

resonant frequency of a short open column or enclosure will give 

1 / 2 

c. Since c = CyRT/M) , the quantity inferred is y, which 
implies a weighted measure of the polyatomicity of the atmosphere. 

In principle, the determination of M and y is made redundant 
by the presence of a mass spectrometer experiment. This redundancy 
is probably worth having for various reasons: (1) Both sensors 

are extremely simple; (2) For a mass spectrometer of limited 
resolution and dynamic range, various isotopes of C, 0, and N 
will complicate the interpretation; (3) Strong inferences can 
still be made about the major atmospheric constituents even if the 
mass spectrometer fails to operate or is deleted from the mission. 

The entire group of sensors should be interrogated at a rate 
which will produce data at least once every kilometer .. .at least 
once in three seconds. Of particular note is the need for 
dynamic range in the pressure and density experiments; P will 
range from about .005 to possibly 100 atmospheres and n from .01 
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19 -3 

to 120 x 10 cm If necessary, multiple sensors with different 

sensitivities and saturation levels must be incorporated. 

Composition of the Atmosphere 

We now consider the major experiment aimed at determining 
the composition of the atmosphere. . .a mass spectrometer. Since 
1954 a number of low-weight non-magnetic instruments have been 
flown on rockets and satellites in the study of the earth 1 s 
upper atmosphere. The difficulties present in the interpretation 
of much of this data have not been minor, and the art may be 
regarded as being possibly in adolescence at this time. Although 
the lower atmosphere of Venus might be thought to be an easier 
experimental subject, by analogy with the Earth’s lower atmosphere, 
no experience exists, to our knowledge, in the mass analysis of 
the troposphere and lower stratosphere from drop-sondes. We 
recommend that the ability of a mass analyzer to determine the 
"exact" atmospheric composition of Venus be regarded very conser- 
vatively until suitable instruments have been built, tested, and 
calibrated under various conditions. 

No specific design recommendations are offered for the 

instrument, Nonmagnetic designs which have been successfully 

flown in upper atmosphere studies include radio frequency, time- 

of-f light, and quadrupole mass spectrometers. Typical resolution 

separates mass numbers up to about 50, and a dynamic range of 

4 

mass current of about 10 is more than adequate to answer the 
first-order questions about the atmosphere. We do suggest that 
mechanical moving parts are quite undesirable and that an entrance 
leak system which takes advantage of the flight profile (into 
progressively high pressures) is to be desired. It is not clear 
that the ability to perform at pressures over a few atmospheres 
would be of any real importance, due to the certainty that the 
lower atmosphere is well-mixed. (Knowledge of composition in 
nonequilibrium microenvironments in the tens of meters 
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immediately above the solid surface is of course desired, but we 
cannot regard this information as being of very high priority at 
the present time.) 


Radiation Environment 

We now consider the class of experiments which measure 
radiation fluxes in different directions and different wavelength 
ranges. Since the number of possible measurements and bits of 
information is very large indeed, selection of a very few 
experiments must be made on the basis of their pertinence to 
answering the most important questions. The particular "mix" of 
sensors suggested here is not in any sense the best that can be 
proposed, but should be typical of the genre. 

Some of the questions which can be answered by photon 
detection are discussed: 

1. Where are the clouds located? - The tops of the clouds 
would be identified by a detector looking directly at the sun in 
the visible. The bottom of the clouds would be detected by com- 
paring the outputs of up- and downr looking detectors (assuming an' 
optically transparent atmosphere below a given level and a surface 
of albedo around 5-30%), The existence of structure, such as well- 
defined cloud decks or cellularity could be inferred in some 
vague sense by having for comparison a third detector looking 
horizontally. 

2 . What is the water vapor profile in the atmosphere? - This 

question should be answered by the mass spectrum. Since we are 
consistently looking for alternative sources of composition data, 
we will consider a spectroscopic method. The most direct method 
is to measure the absorption of the solar flux in a water vapor 
absorption line. The solar flux is strong at 2 microns, but 
scattering by the clouds hurts the interpretation at a shallow 
depth in the clouds . A 6 micron experiment would work deeper in 
the clouds (absorption due to scattering would be down by 3 = 81) . 

At more than one or two optical depths in the clouds, the absorption 
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measured in a band is determined by the mean photon path in 
multiple scattering; in principle, one could make an indirect 
determination of H 2 0 in the clouds by comparing the 2 and 6 micron 
absorption with side-looking and solar-looking detectors. The 
law of diminishing returns has set in, and the mass spectrum is 
the best bet for measuring water vapor in the lower atmosphere. 

The independent determination of water in the upper atmosphere is 
still of real value, since a knowledge of its vapor pressure at 
the cloud tops should tell whether or not the clouds are ice. A 
direct measurement of the cloud composition is quite difficult, 
and redundance of this sort is highly desirable. 

3. What is the radiation balance in the atmosphere? - This 
question is general, and aimed at the theoretical understanding of 
the atmospheric dynamics . The measurement appropriate to the usual 
radiative transfer models is to determine the intensity of the 
upgoing and downgoing radiation streams throughout the visible and 
infrared as a function of wavelength, Disregarding this is 
altogether too grandiose, we suggest experiments in two bands: 

(1) the broad optical band, to detect the total flux in the 
solar black body curve, using an ordinary photodetector in the 
range .35 - 2.0 microns; (2) the broad thermal infrared, to 
detect the radiation streams due to thermal emission by the 
planet and the atmosphere, using an ordinary bolometer, in the 
range 4-20 microns. For each band, an upward looking and a 
downward looking detector with large solid angles of acceptance 
(to include the sun) , and a side looking narrow-angle detector 
(to determine scattered photons) are suggested. 

If all the above suggestions are brought together, we find 
that optical band detectors are relevant to both the cloud 
location and the radiative balance problems. Due to the 
importance of scattering in determining the radiation field within 
clouds, side looking sensors are of general use. A fairly simple 
measurement of absorption in a water band is of real use above the 
clouds, dubious at lower levels. Detection of C0 2 absorption 
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above the clouds in the same way is of course possible. Since 
CO 2 information is implicit in the mass spectrum, the speed of 
sound, the density, and the infrared balance measurements, we 
suggest that any additional redundance is probably not 
worthwhile. 


Composition of the Clouds 

The clouds are postulated to be water, CC^j "dust", and 
various aerosol mixtures which involve ’'dust" and other 
constituents. The lifetimes and falling rates of droplets and 
particles in the size range .1-10 microns are wholly consistent 
with any composition. Considering the very low particle 
densities needed to produce opacity on a several km path, it is 
out of the question to consider a direct sampling. The only 
technique which has come to our attention that would have any real 
chance of working involves detection of cloud particles by counting 
scintillations in a sample tube through which the air is passing. 

By passing the stream into higher temperatures, a second 
scintillation count would determine the amount of evaporation. 

This experiment is specific for clouds. CO 2 clouds can be 

inferred if the temperatures are low enough and the H 2 0 test is 
negative. With a negative H 2 0 test and temperatures above 220 °K, 
dust can be inferred. It is evident that careful temperature 
regulation is needed, and that the experiment is not a very "clean" 
measurement. Since the cloud composition is an important 
first-order question, every effort must be made to devise a 
workable test, of this sort or any other. 

DISCUSSION AND RECOMMENDATIONS 

The atmospheric entry probe and experiment set discussed here 
are derived from a "pure" approach, by considering the major 
questions about Venus and looking for a workable set of crucial 
experiments. From this point of view our major recommendations are 
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(1) The purpose of an entry probe must be the study of the 
lower atmosphere (p > 5 mb) and the clouds. Study of the upper 
atmosphere and ionosphere, while representing a technology 
traditionally connected with the space program, implies the 
measurement of quantities in the Cytherean environment which are 
at best only imperfectly known for the Earth. A real sense of 
priorities demands that the spacecraft and probe be designed as 
far as possible for the prime mission of lower atmosphere study. 

(2) Interpretation of data- from the probe experiments is an 
indirect process, and the different results are strongly inter- 
related, (e.g. mean molecular weight as a clue to major constituents, 
and water vapor at cloud tops as a clue to cloud composition) . 

The integration and intercalibration of the experiments, and the 
matching of these to the needed environment (primarily a thermal 
problem) are a formidable systems design problem. 

(3) For this reason the planning of the probe system and the 
construction and testing of prototype experiments should begin as 
soon as possible. 

(4) Since the probe will be essentially unique, the 
selection of an integrated science package should precede the 
specification of a fixed engineering environment, subject to as 
few constraints as possible (probable ones being weight, general 
shape, and communication capability). 

Our recommendations should be viewed as follows in terms of 
current plans for planetary exploration: 

(1) If a Mariner class of spacecraft can be devised as an 
entry probe, then we urge consideration of a probe of the sort 
discussed here. If a Mariner probe is possible only at the 
expense of many of the proposed experiments, we suggest that the 
minimum configuration which could be tolerated would do the 
following : 

a) Determine the pressure, density, and temperature 

just above the solid surface; 
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b) Measure the atmospheric mass spectrum just above the 
clouds . 

(2) If the Voyager program goes as scheduled, then an early 
entry probe along the lines of this report is strongly urged. We 
recognize that some questions may arise between an early Mariner 
probe or a somewhat later Voyager probe. The resolution of this 
question must depend upon specific times, weights, etc. 
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ABSTRACT 

The object of microwave, millimeter, and infrared spectro- 
scopy has been to provide information about the atmopshere of 
Venus. Measurements at microwaves and millimeter radiometry have 
provided data about the temperature of the surface of Venus of 
about 600°K. In addition, radar measurements have given the best 
value of the dielectric constant of approximately e = 5 +_ 1. 
However, better data can be obtained if further measurements with 
improved equipment are taken in the wavelength region between 3 
and 30 cm. The limit of error can be reduced. The nonresonant 
character of the microwave radiometric data suggests pressure 
broadening of the emission lines. However, infrared data at 
shorter wavelength would show a definite resonance, and in the 
wavelength range of 1 to 10 microns the quantitative nature of 
the C0 2 , CO, H^O, and even other weaker possible components 
could be determined from orbiting satellites about Venus and 
Earth-based instruments on high altitude, airborne or near 
satellite vehicles. 
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Benjamin Lax 
November, 1967 

INTRODUCTION 

Among our neighboring planets the atmosphere of Venus has 
been one of the most intriguing puzzles. To date almost all of 
our information has been obtained from ground-based observations. 
Unlike Mars, which has now been studied by a fly-by spacecraft 
whose occultation measurements have confirmed and sharpened the 
terrestrial data, we still have a great deal of uncertainty con- 
cerning the possible constituents in the atmosphere of Venus. 

Even the magnitude of the total pressure of the components is 
unknown, much less the partial pressures. Among the likely 
candidates, which consists of C0 2 , H 2 0, N 2 , and CO, in addition 
to some small concentration of inert gases, only C0 9 is specified 
with some certainty quantitatively, namely about 10 cm atm 
(NTP) as deduced from spectroscopic observations ^ . It is assumed 
that C0 2 gas is much below the cloud layer, although some 
investigators have proposed that the clouds may have been C0 2 
particles in crystal form. The temperature deduced from radio- 
metric measurements at different regions of the electromagnetic 
spectrum is too high to substantiate such a hypothesis. In regard 
to water vapor there is a great deal of controversy as to its 

presence, although the balloon experiments of Strong and 
2 

co-workers in the near infrared indicate traces of water vapor 
in the Venusian atmosphere. There is less certainty regarding 
the other two likely candidates, namely N 2 and CO. The problem 
in each case, therefore, is to devise techniques within the present 
state of technology to determine (1) the composition of the 
atmosphere, (2) the quantitative features of the constituent gases, 



(3) the structure and distribution of the constituents in the 
atmosphere. In addition to this type of information, radar and 
radiometric data can, in principle, determine the dielectric 
constant and the temperature of the surface of Venus in addition 
to providing qualitative information about the atmosphere of the 
planet. However, even the quantitative nature of the dielectric 
constant is now debated. In analyzing the microwave radiometric 
studies of Mariner II, Pollack and Sagan deduce an absorbing 
cloud of liquid water droplets as the primary source of capacity. 
If the dielectric constant is assumed to be e = 3.6 at 19 mm, 
then they propose an adsorbing and scattering atmosphere contain- 
ing dust , hailstorms or large raindrops with particle radii of 

0.5 mm. The most recent concensus, which is based on the con- 

u 

sistence of the microwave and infrared data , favors water vapor 
and ice clouds as against dust particles as the cause for the 
opacity of the Venusian atmosphere. 

MICROWAVE SPECTROSCOPY 

4 

Microwave measurements have been carried out by radar and 
radiometric techniques from a few centimeters to about 10 meters 
with varying degrees of accuracy. The radiometry has actually 
been performed over a modest band in the 1.35 cm region of the 
water vapor absorption line without definitively indicating 

5 

resonances . Thus, to date the microwave measurements have been 
substantially a method of nonresonant spectroscopy. However, 
these measurements have been very valuable in defining critical 
problems which indicate future directions . The radiometry as a 
function of wavelength has clearly defined two ranges of 
temperature. At the longer wavelength of several centimeters, 
the temperature that is deduced is of the order of 600°K and is 
believed to be that of the surface of Venus. In the region of 
one centimeter or less the temperature drops to below 400°K. 

The interpretation of this is that the shorter wavelengths are 
measuring the temperature of the clouds in the upper atmosphere 
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of Venus. Another- important confirmation is the radar reflection 
which is substantially higher in the longer wavelength region at 
10 cm and above and substantially reduced at 3.8 cm. This is 
believed to indicate atmospheric absorption by gaseous molecules , 
although an alternate mechanism of scattering and absorption by 
particulate material is also a possibility. Much theoretical 
speculation for resonant and nonresonant absorption has been 
stimulated by the former. Unfortunately, the radar data has such 
large limits errors, particularly at long wavelengths, that 
theoretical analysis will not differentiate between the two types 
of absorption much less identify the constituents. However, 
better data may bring about greater consistency with one model 
of the composition than another and reduce the uncertainty of the 
total absorption. The same approach may be applied to the 
determination of dielectric constant which is obtained both from 
radar data and microwave radiometry . The value of the estimates 
have varied from about 2 to 7 with the radar favoring a higher 
value. The uncertainties in the data combined with atmospheric 
absorption may account for the large discrepancies. Consequently, 
it is important to refine and extend these measurements since the 
microwave measurements are the only tools today that can provide 
information about the solid surface of Venus. Now that the tem- 
perature of 600°K has been verified by the successful soft landing 
of the Russian capsule on the surface of the Moon, ground-based, 
airborne, and space-borne microwave instrumentation has an 
excellent opportunity for increasing our qualitative and quanti- 
tative knowledge of the Venusian surface and atmosphere. 

Let us consider the radar reflectivity data as a function of 

wavelength. This is shown in Fig. 1 as reproduced from Evans 
6 

and Pettengill with their indicated error bars. The data shows 
that the 3.8 cm measurements in which the reflectivity is small, 
nevertheless is most accurate. If we then accept this information 
and postulate that at long wavelengths the radar reflectivity is 
determined by dielectric constant of the surface which is assumed 
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to be relatively smooth (the roughness parameter is small) , then 
we can draw horizontal curves at long wavelength. The reflectivity 
is given by: 


r - i£L - 1)2 

" " (/e + l) 2 


( 1 ) 


where e is the dielectric constant of the surface. These are 
indicated on the righthand side of Fig. 1. If we further assume 
that at shorter wavelengths the absorption is due to the tail 
end of a Lorentzian resonance associated with rotational transitions 
of molecular constitutents , the absorption is proportional to the 
square of the frequency. Of the likely constituents the closest 
resonance is that of water at 1.35 cm all other possible 
resonances occur in the millimeter and submillimeter range. 

Hence, the reflectivity can be represented by a simple expression 
of the form: 


a 

R = R m e A (2) 

The ratio R/R„ is most often written as o/o* the ratio of 

oo Q 

the effective radar cross section where the cross section 
2 

Cq = tt a is that of a sphere of radius a with a smooth surface. 

Hence as in the long wavelength limit, in the case of Venus the 

reflection depends on the dielectric constant as in Eq. (1). At 

shorter wavelengths the effective radar cross section or reflectivity 

is reduced by the absorption in accordance with Eq. (2). But we 

know one point with greater certainty than the others; namely, 

R at 3.8 cm. If we take this and the value R . we can draw a 

00 * 

family of curves for different values of the dielectric constant. 
(These curves do not reflect the uncertainty of the 3,8 cm point.) 
From the plot it appears that the dielectric constant of 2,2 from 
the microwave polarization data is outside of the error limits. 

The radar data would indicate from these curves that a dielectric 
constant of 4 to 6 would be a reasonable value or e = 5,0 + 1.0. 
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Figure 1. Radar Reflectivity of Venus. 
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More important, the results indicate that with present techniques 
data could be obtained which would yield such a curve from about 

4 cm to 23 cm with considerably reduced error. This would deter- 
mine both the dielectric constant and the overall nonresonant 
absorption within approximately _+ 10%. 

It is apparent that neither the microwave radiometry nor 
the radar are at present able to detect resonant absorption of 

5 

gases. Barrett and Staelin searched for the 1.35 cm water vapor 
line, but were not able to resolve. This could mean that it was 
not present at the cloud height corresponding to 1 cm emission 
region or else the line was sufficiently broadened to wash out 
the resonance. However, such a radiometric experiment at about 
1.7 mm and near 0.9 mm would show reasonable resonance emission 
if such measurements were attempted from an orbiter, a satellite, 
or airborne radiometer in the millimeter region, and if water 
vapor were present in the predicted concentration from infrared 
data. Such a radiometer may also detect presence of 0 ^ at about 

5 mm and 2 . 5 mm which exhibits sharp rotational lines, although 
the expected concentration in the atmosphere of Venus is small. 

INFRARED SPECTROSCOPY 

It is generally accepted that infrared spectroscopy is a 
more powerful tool than microwave or millimeter spectroscopy for 
identifying molecular constituents of a planetary atmosphere. 

Many of the molecules that have been speculated to be present in 
Venus exhibit strong infrared absorption bands which are quite 
unique for each of the species. Furthermore, in the wavelength 
region of about 1 micron to 16 microns , unlike in the microwave 
region, pressure broadening of these bands will not affect their 
intrinsic character. Of the various compounds that have been 
suggested for Venus, only CO 2 has been positively identified pri- 
marily from its very strong infrared absorption bands at 2.7 
microns and 4 . 3 microns . Among the other candidates which have 
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been proposed, H-0 is one that exhibits many resonant bands. 

L o 

Strong and co-workers have studied the 1.13 micron band and 
7 

Dollfus has observed the 1.4 micron band of H 2 0. In addition, 
water vapor shows even stronger and increasingly broader infrared 
bands centered at 1.8 microns , at about 2 . 7 microns and a very 
broad band from 5 to 7 microns . Fortunately it has a window at 
4 microns and hence would not obscure the unique band of C0 2 in 
this region. Another candidate which naturally may accompany the 
presence of CC> 2 is CO. This has a strong identifying band at 
about 4.7 microns which should be detectable, although the anti- 
cipated concentration suggests a weaker absorption. However, with 
an orbiter or improved ground-based instrumentation this should 
now be possible. Nitrogen itself, which has been considered, 
does not possess either a microwave or an infrared spectrum, but 
its compounds N 2 0 and N0 2 have both, N 2 0 can be distinguished 
from CO in that, in addition to a band at 4.7 microns, it has one 
at 7.8 microns on either side of the large water vapor band, 
making its spectrum detectable. 0 2 and 0^ both have infrared and 
microwave resonance spectra. 0 2 has weak transitions at 1.06 and 
1.27 microns and stronger ones at 0.762 and 0.688 microns. Lastly 
among the molecules with strong infrared spectra, CH^, one of 
the likely hydrocarbons , does exhibit relatively narrow bands at 
3.3 and 7.7 microns again in a fortunate position relative to 
water vapor. 

Recent high resolution measurements by Connes , Benedict, and 

g 

Kaplan , using a Michelson interferometer, has shown unique 
spectra of Venus. These were taken at the Observatorie de Saint 
Michel at high altitudes. The spectra were taken in the three 
atmospheric windows of water vapor between 1 and 3 microns. The 
data exhibited distinct lines attributed to HC1 and HF, and even 
distinguished between isotropic species of Cl. This particular 
result demonstrates the extreme value of high resolution spectro- 
scopy in the infrared. Such an instrument above the Earth’s 
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atmosphere or in a Venus fly-by or preferably an orbiter would 
uncover other components of the upper atmosphere of Venus. 

GENERAL COMMENTS 

Among the important major components, N£ is one of the most 
difficult to detect because it has no active resonance spectrum 
either in the infrared or the microwave-millimeter regions 
because of its symmetrical charge distribution. The inert gases 
such as argon and neon which may be present are even harder to 
detect. The only likely possibility will be to look for their 
natural lines in the ultraviolet. This is a more difficult task 
and in the presence of either molecules which also absorb in the 
ultraviolet this becomes even more involved. 

In examining Fig. 2, which shows pictorially some of the 
spectral vibrational-rotational bands of interest, some pertinent 
observations can be made. The ground-based observations, which 
have been made with the Mount Palomar telescope and Mariner II 
for the observation of Venus in the infrared, were made essentially 
in the 8 to 14 micron window of ^0. With the exception of 
this does not include spectral lines of importance to us, The 
negative results may indicate that 0^ is not present at all or not 
in sufficient quantities to be observed radiometrically. Hence 
it appears that it would be of greater significance to do measure- 
ments in the spectral range of approximately 1 to 10 microns which 
include the water vapor lines. However, it is necessary to get 
sufficiently high above the earth's atmosphere to eliminate its 
water vapor. A limited advantage may be obtained by observations 
on mountain tops in dry regions, such as Arizona. However, much 
more can be achieved by the use of high altitude balloons or 
satellites. The possibility of using high flying airplanes with 
appropriate large area mirrors is another suitable alternative. 
Finally, an orbiting radiometric scanning spectrometer about 
Venus in this region of the infrared would be most desirable. 
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Absorption 



Figure 2. The near-infrared solar spectrum (bottom curve). 

Other curves are laboratory spectra of the molecule, 
indicated. The diagram is from J. N. Howard, 

D. L. Burch, and D. Williams, Geophysical Research 
Papers (U.S.) 40_, Rept . AFCRL-TR-5 5-213 (1955). 
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The priority should be in terms of relative cost which would 
probably be in the order of planes or balloons (or both), then 
satellites orbiting about the Earth, and finally the Venusian 
orbiter or fly-by if it precedes an entering capsule which 
successfully samples the planet's atmosphere. 

Greater detail is precluded from this report since its pur- 
pose is not to spell out the specific instrumentation, but to 
point out the fruitful regions of microwave and infrared spectro- 
scopy which can be exploited within the next five or ten years . 
The recommendations are based upon improved capabilities where 
some success already has been achieved as in the microwave 
region. In the infrared region, technology of detection and 
spectroscopy in the 1 to 10 micron region gives added incentive 
to try airborne and spaceborne measurements . 


Added Note 

This report was originally written prior to the successful 
launching of the Soviet Venus landing probe. This established 
that CC >2 was the major constituent of Venus and that the surface 
temperature was very high. However, the high temperature pre- 
vented the electronic equipment from operating for a sufficiently 
long period to give all the information that is desired. Hence 
this would still justify future studies by orbiting satellites 
about Venus and Earth-based instrumentation whether on the ground 
or above the Earth's atmosphere. 
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SOME COMMENTS ON THE SCATTERING OF THERMAL MICROWAVE 
RADIATION BENEATH A PLANETARY SURFACE 


ABSTRACT 

The microwave radiation thermally emitted by a planet 
originates in the planetary medium and reaches the surface by a 
process of radiative transfer. A simple model is used to 
determine the effect on the transfer process of single scattering 
from centers randomly distributed in the medium. Particular 
reference is made to the ability of the scattering mechanism to 
impose a weak linear polarization upon the radiation. It is found 
that the polarization from scattering is much weaker than that 
induced by transmission through the surface. As it stands the 
model is incapable of providing an explanation of the low value of 
the surface dielectric constant of the Moon and Venus deduced from 
measurements of the microwave polarization. The suggestion is 
made that the development of a multiple scattering model is 
necessary to provide a complete understanding of the interaction 
of microwave radiation with a planetary surface. 
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SOME COMMENTS ON THE SCATTERING OF THERMAL MICROWAVE 
RADIATION BENEATH A PLANETARY SURFACE 


B . G . Smith 
July, 1967 

INTRODUCTION 

The Moon and planets emit microwave radiation and an analysis 
of the characteristics of the radiation can give information about 
the physical parameters of the source. The microwave emission 
from the Moon, Mars and Venus is of low intensity, with weak 
linear polarization and, apart from solar phase effects, time 
invariant. The simplest explanation is that the radiation is from 
an incandescent black body. 

The model to describe thermal microwave emission from a 
planetary surface was first elaborated by Troitskii (ref. 1). 

(For convenience in this paper the Moon will be included under 
the general heading of Planets, whenever no confusion will arise). 
Radiation leaving the surface originates in the planetary medium 
up to an absorption length deep. The intensity of the radiation 
is a measure of the brightness temperature of the medium averaged 
in some way over this depth. Microwave radiometers in present 
use on Earth have insufficient angular resolution to allow 
examination of local features on the planets and in fact sample 
the radiation integrated over the whole planetary disc. Average 
disc brightness temperatures have been derived for the Moon, Mars 
and Venus. Those for the Moon and Mars are compatible with the 
measured infrared temperatures, but a well known discrepancy 
exists between the two temperatures of Venus, providing fruitful 
ground for atmospheric model makers. 

A second characteristic of thermal radiation is its polar- 
ization. If the planetary medium is isotropic and if scattering 
plays no part in the transfer process by which radiation reaches 



the surface, the upward stream of radiation beneath the surface 
will be unpolarized, Passage through the surface will then impose 
a linear polarization upon the wave, If the surface may be 
considered locally as a plane interface between two dielectric 
media , the magnitude of the induced polarization will depend upon 
the angles of incidence and refraction, according to Fresnel's 
transmission equations. The total radiation integrated over a 
uniformly emitting spherical surface would, however, from 
arguments of symmetry, be unpolarized. Since the sun is the 
ultimate source of at least part of the radiated energy , the 
variation in solar phase angle produces a surface temperature 
variation with planetary longitude and this in turn causes 
different parts of the surface to contribute different amounts 
to the total radiation. The integrated emission will in this 
case be weakly polarized. In addition, the weighting accorded 
different contributions depends upon the characteristics of the 
radiometer used to make the experimental measurement. The 
orientation of an interferometer axis defines a preferred 
direction of the planetary disc preferentially weighting the 
contributions from the surface lying on a diameter perpendicular 
to the axis. In this case, too, the integrated emission is weakly 
polarized. A twin-antenna microwave interferometer has been used 
in this way to study radiation from Venus (ref. 4). 

Fresnel's transmission equations contain the dielectric 
constant of the subsurface medium as a parameter. Measurements of 
polarization have given values for the dielectric constant, e^, 
of the surface layer of the Moon and Venus quoted in Table 1. 

There they are compared with values of deduced from the radar 
backscatter cross section of those bodies. The discrepancy is 
striking, and indicates that in each case the polarization of 
thermally emitted radiation is less than that expected from a 
surface layer with the radar dielectric constant. Various 
explanations have been offered. Roughness of the planetary 
surface would impose less polarization upon the emergent 
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radiation than the equivalent smooth surface. Hagfors and 
Moriello (ref. 8) have examined a rough surface model and are 
able to explain some, but not all, of the discrepancy. Another 
explanation is that the planetary medium is inhomogenous with 
depth, with compacted material of higher dielectric constant 
underlying a more diffuse surface layer (ref. 5), It is 
suggested that radar will sample the dielectric constant of the 
lower strata while the polarization of the thermally emitted 
radiation will be induced by passage across the outer boundary 
of the diffuse top layer. 

This paper describes a new model which incorporates 
scattering of the radiation by subsurface inhomogeneities in 
the process of radiative transfer through the planetary medium. 
Scattering mechanisms themselves can impose polarization upon 
initially unpolarized radiation. Depending upon whether the 
polarization from scattering reinforces or subtracts from the 
polarization induced by transmission through the surface, the 
dielectric constant, according to the simple model, will be 
over or under-estimated. The possibility of explaining the 
dielectric discrepancy in terms of a scattering model provided 
the stimulation for the present study. In anticipation of the 
conclusion it must be stated that in fact the effect of scattering 
is to produce a polarization of incorrect sign and of insufficient 
magnitude to account for the discrepancy. 

DESCRIPTION OF THE MODEL 


Scattering will play an important role in radiative transfer 
if the scattering length Z of radiation in the medium is 
comparable with or less than the absorption length Z^. A crude 
model will be used to estimate the effect of scattering centers 
embedded in the medium. The planetary surface, illustrated in 
Fig. 1, is locally plane, separating the planetary medium of 
dielectric constant from the . atmosphere (if any) of dielectric 
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flux I ( u ,z) 


Figure 1. - A Section of a Planetary Surface 


constant (id). The medium is homogeneous except for the 
presence of a random distribution of inclusions of dielectric 
constant C/e^), of characteristic dimension less than a free 
space wavelength X. The interaction of radiation with an inclusion 


1-5 


will be approximated by Rayleigh scattering, for which the total 

scattering cross section a is (ref. .9.3 

s 



where a is the polarizability of the inclusion. For spherical 
scattering centers of radius R, a is given by (ref. 10) 


- £_ 


a 


2 £, + e. 


R 3 e. 


( 2 ) 


The scattering length £. expressed in terms of the scattering cross 

section and the number density of scattering centers n , is 

s 


i 

s 


n a 
s s 


(3) 


4 7T 3 

If g ( =— R n g ) is the fractional volume of the medium occupied 
by spherical inclusions 



It is not difficult to imagine an inhomogeneous medium for which 
typical values of the parameters would be 6^0.05, R^O . 1A , e.^4, 
£^2 and in this case Jl^lOA. Estimates of the absorption length 
of microwave radiation beneath the lunar surface have been made 
(ref. 11); a typical value is 1 10A . Although criticism can be 

raised about the validity of this estimate, the deduction may be 
drawn that there is no strong a priori reason why scattering of 
radiation should not be an important process within the lunar 
medium. It will be assumed that a similar argument could be 
advanced for the planets. 



The equation of radiative transfer within a plane parallel 
planetary medium (ref, 9) is, referring to Fig. 1, 


„ dl (y , z ) 
v dz 



I(U,z) 



B(z) 


+ 


1 

l 




t t » 

P(y<|> s y <P ) l(y jz) 


(5) 


I (y,z) is the radiative energy flux per unit solid angle at a 
depth z in the medium in a direction (9,0) making an angle 0 
with the surface normal. y = cos 0, and because the system has 
cylindrical symmetry about the z axis I is independent of azi- 
muthal angle <j> . The term on the left of the equation represents 
the change of the energy flux in traversing a path length (-dz/y) 
in the medium (z being measured downwards). On the right, the 
first term represents attenuation of the flux by scattering and 
absorption. The second term describes. the isotropic emission 
process according to Kirchoff's Law. B(z) is the Planck 
function at the appropriate wavelength and depends upon temperature 
which may be a function of depth. The third term represents the 
enhancement of the flux caused by scattering into unit solid angle 

t » 

in the direction (0,0), and P(y0 ,y <p ) is the angle dependent part 
of the angular scattering cross. section. Since the radiation can 
be polarized, it is necessary to differentiate between a flux 1^ 
polarized with its E vector in the plane containing the direction 
of the flux and the. normal to the surface , and its orthogonal 

. /M 

counterpart • I is thus a two component quantity I ^ J, and P 


a two by two matrix, which for Rayleigh scattering has the form 
(ref. 9) 
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where Q is a matrix with the property that 



Q<8MV) ~ 0 


(7) 


Q plays no part in the following analysis. B is also a two- 
component quantity, but for an isotropic medium the components 
and B2 are identical. 

If scattering is a less important attenuative process than 

absorption, equation 5 can be solved by a series expansion in 

powers of . The first two terms in such an expansion 

describe the transfer of radiation suffering zero and one 

scattering process. In order to demonstrate the physical 

properties of the model more clearly higher order scattering 

processes will be neglected. Developing the series from an 

intuitive point of view, there will be two constituents of the 

radiation reaching the surface, illustrated in Fig. 2. Radiation 
(A ) 

I emitted by each volume element reaches the surface unde- 
flected after suffering attenuation by absorption and single 

(A) 

scattering en route. I satisfies the reduced equation of 
transfer 
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radiation suffering attenuation by 
absorption and scattering while in transit 

radiation suffering attenuation by 
absorption only 


Figure 2. - Constituents of the Radiation Reaching the Surface 
for the Single Scattering Model 
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Equation 8 may be integrated to give 



-U/* A * U '- z) 

B(z')e 

-<1/4. + l/4 = )i (z'~z) 
A s y 

B ( z ) e 

(9) 


where the boundary condition has been imposed that there be no 
downward radiation stream at the surface (z = 0), 

( B) 

The second constituent of the radiation, I , reaches the 
surface after suffering a single scattering, and attenuation by 
absorption before and after the scattering process. I satisfies 
an equation of transfer 


-U 


dl 


(B) 


dz 


1 T (B) , . 

■T 1 + * 


( 10 ) 


( B ) 

where $ ( y, z) , the two component source function for I , is 


$ ( u , z ) 




» I a 1 » 

p(e<t>,e ^ ,z)(ii) 


—(A) (A) . . . 

and I is I evaluated xn the lxmxt SL -*■ 00 , since the source 

( B ) ® 

of I is radiation which has suffered no earlier scattering. 

Implicit in the definition of $ is a neglect of radiation which 

suffers single scattering after a reflection from the surface. Its 
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contribution is small and for ease of calculation will be excluded. 
Equation 10 can be integrated to give 1^®^ 

(B) 1 f -z/u^A 

U>0; I (y ,0) = — I dz f $(y,z')e A (12) 

If the brightness distribution with depth, B(z), is specified the 
integrals in Equations 9, 11 and 12 may be evaluated to give I (A) 
(y,0) and 1^^ (y,0). Their sum, l(y,0), for y>0, is the upward 
flux of radiation at the surface within the medium, which sub- 
sequently suffers refraction and a transmission loss on crossing 
the surface. If J k (y) is the two component radiative energy 
flux per unit solid angle outside the surface in a direction 
(0,<j)), making an angle e" with the surface normal, and y = cos o' 

y>0; J k (y) = (1 - R k ) I k (y,0) k = 1,2(13) 

where k is an index denoting one of two polarization components 
and sin 0 = /e^ sin 9 according to Snell’s law of refraction. 

R k is the appropriate Fresnel surface reflection coefficient 
given by (ref. 12) 


R 


1 


tan 2 (9-6 ) 
tan 2 (6+"0 ) 


R, 


sin 2 (8-6 ) 
sin 2 (0+0" ) 


(14) 


The factor l/e 1 appears in order to allow for the change in solid 
angle occurring at refraction through the surface (ref. 1). 

The purpose of this paper is to compare the polarization 
induced by scattering with that induced by transmission through 
the surface. It is convenient to define the polarization y 

s 

caused by scattering alone. 
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(15) 


I 1 (U,0) - l 2 (y,0) 
n s = I^y/ti) + l 2 (li,0) 


y>0 


and to compare it directly with the polarization caused by 
transmission through the surface in the absence of scattering 

J 1 <v) - J 2 (y) 

n T = z 

J^Cy) + J 2 (y) 




(16) 


In this limit, I, (y,0) = (y,0), y^.0 , and from equation 13, 

(1 - R.) - (1 - R-) 

n T = — 

(1 - R x ) + (1 - R 2 ) (17) 


Explicit evaluation of n will be made for two brightness 

S 

distribution functions B(z): that of an isothermal medium, and 

a linear brightness gradient. 

(1) Isothermal medium . B(z) = constant, B^sB^ = Bq 

(A) 

Equation 9 can be integrated directly to give I . At the 
surface 


y>0 ; 


T (A) 
x k 


(y ,0) 


Vl 

‘V*Aj 


B, 


k = 1,2 (18) 


a result which, in the limit of no scattering, reproduces that of 
Troitskii (ref. 1). Using equations 9, 11 and 12 the expression for 
1^ - I 2 in the isothermal medium is 


y>0 ; 


l 1 (y,0) - l 2 (y,0) 


2£ 1 B 0 



f 


(y) 


(19) 
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where 


r 




f(u) = 


P (1-p ) 


(1-3 p; log (1+1/ p) - (i- - 3 p) 


( 20 ) 


The total intensity in the two polarizations is little affected 
by scattering, and is given approximately by 


P>0; l 1 (p,0) + l 2 (p,0) sr 2e 1 B Q 

with the result , from equation 15 , 

JU 

n s ~— f(y) 

*s 


( 21 ) 


( 22 ) 


n g (for ~ 1) and r^, calculated for two values of are 

displayed as functions of "e in Fig. 3. Two features stand out. 
One, n s has the same sign as a nd so for a real surface the two 
polarization mechanisms will reinforce one another. Two, n is 
at best an order of magnitude smaller than rVp. 

1 

(2) Linear brightness gradient . B^(z) = B 2 ( z) = B^+B z 

With this expression for B(z) the integrals in equations 9, 

11 and 12 lead to 


p>0; I, (p,0) - I 0 ( p , 0 ) = 


2e l It 1 ) CB n + B ^ f |^)(1-P 4 ) B A a e 


32 A. 


A 1 


(23) 


A simple interpretation exists for equation 23. The first term 
to the right of the equality sign is equivalent to the contribution 
from an isothermal medium of brightness equal to that at a slant 
depth of one absorption length. The second term, which is of 
opposite sign, is equal to the difference in the intensities of 
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Figure 3. -A Comparison of n and Tim as a Function of 0 

S X 

Evaluated for JI./& = 1, e, = 2 and 4 (Remembering 

S X 

\i = cos 9 j y = cos 9 and /IT sin 9 = sin 9)* 




the two polarized components of scattered radiation produced when 
a steady vertical flux of total intensity j B interacts with 

the scattering centers contained in one absorption length depth of 
the planetary medium. The second term is important only if the 
brightness gradient is steep, in particular only if the brightness 
more than doubles in the distance of one absorption length below 
the surface. This is unlikely to be the case for a planet, and 
so the second model is essentially . indistinguishable from the 
first . 


CONCLUSIONS 

Disappointingly, it is impossible to understand the dielectric 
discrepancy in terms of a simple scattering model. Indeed, if such 
a model were applicable as it stands, current estimates of from 
polarization measurements should be slightly reduced rather than 
increased. 

The calculation does, however, provoke thought about the 
relative importance of absorption and scattering as attenuative 
mechanisms beneath a planetary surface. The current estimated 
value of the absorption length for microwave radiation in the 
lunar medium (ref. 11) is based upon an interpretation of the 
emission process in terms of radiative transfer without 
scattering, an interpretation which this paper is calling in doubt. 
If the absorption length is in reality considerably greater than 
this estimate, it will be necessary to perform a completely new 
analysis of a multiple scattering model for the process of 
transfer by which microwave radiation diffuses through a planetary 
medium. An analysis of this sort might have interesting conclusions 
for the interpretation of both microwave emission and radar 
backscatter experiments from planetary surfaces. 
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WATER ON MARS 


ABSTRACT 

Water in the form of ice can exist on Mars as permafrost 
which is either in equilibrium with the water content of the 
atmosphere or is gradually evaporating through a protective layer 
of soil. The latter situation is evaluated quantitatively for 
grains from 0.5 to 1000 microns and for porosities from 0.01 to 
0.8. The effective diffusion coefficients are calculated and the 
required thicknesses of the protective layers estimated. It is 
concluded that such subsurface ice can exist at Martian latitudes 
lower than those for which the permafrost is in equilibrium with 
the atmosphere. This may lead to seasonal variations of radar 
reflectivity and may explain the higher reflectivity of the dark 
areas as compared to the bright areas. 
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WATER ON MARS 


R. Smoluchowski 
July, 1967 


INTRODUCTION 

No doubt there is water on Mars ^ . The biological, 

climatological and morphological implications of this fact depend 

radically upon its actual amount, form and location. The only 

direct quantitative measurement available is that of precipitable 

water vapor in the atmosphere in the amount of 14 + 7 microns. 

( 2 ) 

The generally subfreezing surface temperatures imply also that 

( 3 ) 

there may be considerable amounts of ice as subsurface permafrost 

( 3 ) 

This has been discussed in detail by Leighton and Murray . The 

polar caps now appear to be solid CC^ rather than ice as previously 

assumed. One can be sure also that there are no open bodies of 

( 4 ) 

water but the occurrence of surface temperatures as high as 

340°K suggests that there may be localized temporary and periodic 

( 5 ) 

melting. In fact the best values of the triple point of water 
6.11 mb and 273°K fall right in the middle of the range of pressures 
deduced from the Mariner IV occultation experiments v and also are 
bracketed by the calculated and observed surface temperatures. The 
occurrence of pressures and temperatures appropriate for liquid 
water, however, is rarer than those for ice and vapor. 

The purpose of this note is to look in a quantitative manner 
into the stability of permafrost, its possible melting and the 
preservation of moisture in the Martian soil. A possible connection 
between these phenomena and the behavior of dark areas and their 
radar reflectivity is also discussed and suitable observations 
suggested. 



EFFECT OF POROSITY ON DIFFUSION 


If the mean annual temperature is low enough, permafrost can 
exist in equilibrium with the water content of the atmosphere. 

Under this condition permafrost on Mars can exist within a few 

( 3 ) 

meters from the surface at latitudes higher than about 45 degrees 

In this paper we are concerned rather with a non-equilibrium 

situation when the mean temperatures are higher and the ice has 

a life time determined by the thickness of a protective layer. 

To a certain extent this resembles the problem of the retention 

of ice on the moon. It is thus important to estimate the rate of 

permeation of water vapor through layers of particulate matter 

which one would expect to encounter on Mars. It is rather 
C 4 7 ) 

certain 9 that the well-known yellow clouds require particles 

not less than a couple of hundred microns in size but these 

particles constitute probably only a small fraction of the actual 

C8 ) 

surface coverage. A careful analysis of the thermal properties 
of Martian surface suggests that the predominant size is from a 
few tenths of a micron up to about 20 microns. The larger 
particles are undoubtedly mixed in but the effective pore size 
will be determined by the smaller grains. 

The effective pore size or the hydraulic radius m depends 
not only on the grain size but also on porosity. We shall con- 
sider here grains with diameter $ equal 0.5, 10, 200 and 1000 
microns and porosites .01, 0.1, 0.2, 0.5 and 0.8. In Table 1, 
standard permeability formulae give the corresponding values of m 
expressed in p : 
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Table 1 


EFFECTIVE PORE SIZE OR HYDRAULIC RADIUS m 


Grain 

Diameter 

Porosity 

e= . 01 

6 = 0.1 

e= 0.2 

e = 0.5 

CO 

II 

0 

* 

CO 

4> = 0 . 5 p 

8xl0 -4 

.01 

0.021 

0.083 

0.332 

<j) = 10]i 

1 . 7xl0 -2 

.02 

0.42 

1.7 

6.8 

<j) = 200y 

0.33 

3.9 

8.3 

33 

130 

<p = 1000y 

. 

1.7 

20 

42 

169 

667 


( 6 ) 

Assuming that the Martian pressure near the surface is 

5mb one obtains in the highly diluted QO^-B.^0 mixture for the 

mean free path X of CO 2 molecule 4-8p for an H 2 O molecule for 

temperatures varying from 160 to 300°K. These values fall in the 

middle of the range of pore diameters listed in Table 1 which 

indicates that the mechanism of transport of water vapor will vary 

from normal diffusion in large pores to Knudsen flow in small 

pores. Under these conditions the very complete theory of 

diffusion in porous media at uniform pressure as given rather 

( 9 ) 

recently by Evans, Watson and Mason has to be used. In the 
normal diffusion region in a binary gas mixture the flux of 
molecules A along z axis is given by 


J A = D AB ( 3 V 3z) + X A J 


( 1 ) 


where D^g is the usual interdiffusion coefficient, n^ is the 
number of molecules A per cc. , is their mole fraction and 

J = J A + J B . 
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In the Knudsen region 


■ D AK <9n A m) 


( 2 ) 


applies where is the Knudsen diffusion coefficient or Knudsen 

permeability for molecules A. In the intermediate region the 
effective diffusion coefficient is given by the Bosanquet 
relationship : 


^ D A^eff = (D AB ) eff + (D AK ) 


(3) 


where 


(D An ) 


V M B 


and 


AB eff 16 nq |M A M fi 
1/2 


2 tt RT 


1/2 




n _ 0 m 2 ttRT 
AK " q "ITT" 


8 +nb 


(4) 


(5) 


In these formulae e is porosity, n is the total molecular 
concentration, q is the tortuosity factor, and Mg are molecular 
weights, S is the scattering cross section, Q is the collision 
integral usually close to unity and b is related to the ratio of 

diffuse to specular scattering and is also close to unity. For 

°2 ( i n 1 

the f^O-CC^ system S = 52 A and also v ' 

D = D (T/T ) 2 (p /p) 
o o -o r 

where we chose T = 273°K and p =1 atm. From Eqs. (4) and (5) 

o o 

one obtains: 


and 


^ D AB^eff = 20,3 q ( T )2 

^ o 


n _ „ Q , n 4 m ,T . 1/2 

d ak = 4 - 8 * 10 5 ( t^ 
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The flux of water molecules is then given by 



(D A } eff 


(9n A /3z) 


where is given by Eq. (3). The second term in Eq. (1) 

has been omitted because of the very small values of X^. The 
only quantity in these formulae for which there are no reliable 
theoretical values is q, the tortuosity factor. This factor is 
supposed to take care of the fact that at low porosities many of 
the pore paths lead to side-wise or even backward diffusion and 
it is usually determined from experimental data. For porosities 
in the range 0.2-0. 5 parameter q turns out to be around 5, for 
lower porosities it is higher. 

It appears that between 160 to 300°K the influence of 

temperature on the effective diffusion coefficient is negligible 

as compared to the influence of the mean pore diameter. Table 2 

2 

shows the resulting mean effective diffusion coefficients in cm 
sec 1 for the various structures of the porous material. It is 
assumed that the tortuosity factor q = 1 for e = 0.8, q = 5 for 
e = 0.5 and q = 10 for the lower porosities. 


Table 2 

MEAN EFFECTIVE DIFFUSION COEFFICIENTS IN CM 2 SEC 


Grain 

Diameter 

Poi 

?osity 

e=. 01 

£ = 0.1 

£-0 . 2 

£=0.5 

£ = 0.8 

<J> = 0 . 5p 

. 0004 

.004 

.008 

. 07 

1.3 

<j> = 10U 

.004 • 

.06 

0.1 

.7 

9 

cf> = 200y 

.009 

.08 

0.2 

1.1 

11 

<|> = 1000y 

.01 

.1 

0.25 

1.2 

12 
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THICKNESS OF THE PROTECTIVE LAYER 


Our model consists of ice under a layer of thickness L of 

porous soil and it is assumed that in order to have appreciable 

permafrost or subsoil ice on Mars a layer of ice 10m thick should 

have existed over a period of a billion years. This leads to a 

-13 2 

maximum loss rate of 3x10 gr/cm /sec. The partial vapor 
pressure of water above ice^^ is given by its temperature T and 
we assume that at the surface, temperature T , the partial 
pressure of water is zero. Thus one has a system in which the 
water vapor diffuses down a concentration gradient and along a 
temperature gradient. In view of the fact, mentioned above, that 
temperature has only a secondary effect on the diffusion coefficient 
this influence will be neglected. 

With these numerical data one can calculate the thickness L 

of a layer of a given porosity and grain size which would be 

sufficient to keep the loss of water vapor from ice at temperature 

- 1 3 - 2 - 1 

T below the value 3x10 gr cm sec . The results are given 

in Table 3. They have primarily relative significance because 

the absolute values are uncertain at least within a factor of 2 

or 3. For e = 0.8 the calculated thicknesses L were greater 

than 10 meters. It should be stressed that the dependence 

of L on T reflects the variation of the pressure of water vapor 

over ice and not the variation of the diffusion rate which as 

mentioned above is only a slow function of temperature. The 

range of temperatures in Table 3 covers the range of mean 

( 3 ) 

annual temperatures on Mars. 

It appears that there is a considerable range of particle 
sizes and porosites for which a layer of ice could be preserved 
during a time of the order of a billion years under a layer of 
reasonable thickness. This implies that ice or permafrost may 
exist at latitudes lower than 45 degrees. 
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Table 3 


CALCULATED THICKNESS L 


Grain 
Diameter 
(<t> in y) 

Temperature 
(T in °K) 

Porosity 

£=.01 

E — . 1 

E = . 2 

e= . 5 

0.5 

240 

10m 


■■ 



230 

3 . 3m 


II 



2 20 

lm 

10 m 

1 



210 

30cm 

3m 

1 



200 

6 cm 

6 0 cm 


9m 


190 

1. 5 cm 

15 cm 

30 cm 

2.2m 


170 


2 . 5cm 

5 cm 

36cm 


160 




6 cm 

10 

220 

10m 





210 

3 m 





200 

6 2 cm 

9m 




19 0 

15 cm 

2. 2m 

3 ... 7m 



180 

2 , 5 cm 

36 cm 


2.5m 


170 


6 cm 

9 cm 

3 6cm 


16 0 



1cm 

4cm 

200 

210 

5 . 8 m 





200 

1 . 2 m 

9m 




190 

30 cm 

2 . 2m 

3 . 7m 



18 0 

4 . 8 cm 

3 6 cm 

6 0 cm 



170 


5 . 4 cm 

9 cm 



160 



2 cm 


1000 

210 

7 . 3 m 





200 

1. 5 m 





190 

37 cm 

3 . 7m 

9 . 3 m 



18 0 

6 cm 

60 cm 

1 . 5m 



17 0 


9 cm 

22cm 



160 


1cm 

3cm 
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SEASONAL AND DIURNAL VARIATIONS OF SURFACE TEMPERATURE 


In discussing the existence of ice and water on Mars one 
has to take into account the seasonal and diurnal variations of 
the surface temperature at various latitudes. At a sufficient 
depth below the surface these variations are negligible but they 
may reach 100° and more at the surface . It follows that the 
temperature gradient between the layer of ice and the surface is 
sometimes positive and sometimes negative. In the first case the 
water vapor will escape as discussed above. In the second case, 
however, it may condense in the colder outside layers and form ice 
which will exist until the temperature gradient is gradually 
reversed. In such pores ice may be very stable since only a small 
fraction of its surface will be exposed to the atmosphere. In this 
connection it should be kept in mind that the results shown in the 
tables above were obtained for a mean pore diameter m given by <p 
and e. There are, however, always some larger and some smaller 
pores. For instance^ 1 ^ for a fixed porosity 0.5 and grain size 
decreasing from about 4mm to less than 0.5mm the ratio of capillary 
volume to non-capillary volume increases from 0.25 to 15. Thus 
even for a large non-capillary value of m there is a considerable 
fraction of capillary pores . Since in the warmer regions of Mars , 
temperatures well above the freezing point were observed, it follows 
that under favorable circumstances some water may be retained at 
least for a while in capillary pores. This may have biological 
consequences . 

The existence of ice in intergranular pores implies that ice 
may also form in cracks of the grains themselves . Such condition 
leads to a gradual break up of the grains and to a comminution of 
the surface material in analogy to well known geological phenomena 
on Earth. 
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EFFECT OF ICE OR WATER CONTENT ON RADAR REFLECTIVITY 

It is interesting to speculate whether ice of water content 

could affect the radar reflectivity of the soil. An analysis of 
• ( 12 ) 

this reflectivity indicates that the effective dielectric 

constant of the bright areas is about 2.2 while for the dark 

areas one obtains 3.5. Typical values for sand are near 2.8 as 

compared to the dielectric constant of water about 80 and of ice 

3.3 to 130 (depending on frequency and temperature). It appears 

( 13 ) 

thus quite possible that the differences in radar reflectivity 
could be interpreted as due, at least in part, to ice or water 
retained in the subsurface layers. The thicknesses of the pro- 
tective layers listed in Table III and those calculated from 

( 3 ) 

equilibrium considerations are quite comparable to the penetration 

depth of the 3.8 to 70cm radar wavelengths used in various Martian 

studies. Without further knowledge of the gradients of porosity 

and of particle size below the surface it is not possible to 

calculate the actual distribution of ice in the soil. Nevertheless, 

even if the radar does not reach the permafrost itself it will easily 

reach the seasonal condensation layer discussed above. According 

to various studies of the thermal properties of the Martian sur- 
( 38 ) 

face 5 , the amplitude of seasonal temperature variations 

decreases by a factor of ten within a layer of the order of several 

centimeters. If this is the case then the shorter wavelength 

reflectivity would be more likely to show seasonal variation than 

( 13 ) 

the longer wavelengths. Recent radar results suggest that 

Martian dark areas are probably lowlands rather than highlands 
and thus one is tempted to ascribe to them higher humidity content 
and a higher reflectivity than to the dry bright uplands. One way 
to check this conclusion would be to measure seasonal and diurnal 
radar reflectivity changes of bright and dark areas. This may be 
difficult to accomplish from the Earth because of the small 
variation of the phase angle and the large variation of distance 
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between Earth and Mars. On the other hand, an orbit er would be 
ideally suited for this purpose. 

Planetary radar reflectivity is a function of many parameters 

and thus its seasonal variation, if it exists, could be accounted 

for on several models. Nevertheless, the ice-water mechanism 

seems more realistic and more attractive than an analogous 

C 14 ) 

mechanism based on variable photoconductivity of the soil. 
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